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Chapter 1 



Introduction 



The Standard Model of elementary particle physics provides a most successful descrip- 
tion of the electroweak and strong interactions among all presently observed particles. 
However, the observational fact that most of the matter of the Universe resides in the 
form of cold non-baryonic dark matter provides an impressive evidence for physics be- 
yond the Standard Model [1]. In fact, the nature and identity of dark matter is one of 
the most pressing questions in the natural sciences. 

Remarkably, supersymmetry (SUSY) offers an attractive solution of the dark matter 
problem. Supersymmetry assigns to each particle a superpartner whose spin differs 
by 1/2. One can distinguish between "normal" matter, i.e.. Standard Model fields, 
and their super symmetric partners by assigning to each particle the ii-parity quantum 



number, and the spin of the corresponding particle, respectively. For conserved i?-parity, 
the lightest supersymmetric particle (LSP) is stable. Thus, the LSP is a compelling 
candidate for dark matter, provided that it does not have electromagnetic or strong 
interactions. 

The gravitino G is a particularly attractive candidate for such an LSP. Any supersym- 
metric theory containing gravity predicts the existence of the gravitino, a spin-3/2 par- 
ticle which aquires a mass from the spontaneous breakdown of supersymmetry. As the 
superpartner of the graviton, it is extremely weakly interacting. Hence, if the gravitino 
is the LSP, it can be dark matter. 

At high temperatures, gravitinos are generated in inelastic scattering processes with par- 
ticles that are in thermal equilibrium with the hot primordial SUSY plasma. Assuming 
that inflation governed the earliest moments of the Universe, any initial population of 
gravitinos must be diluted away by the exponential expansion during the slow-roll phase. 




where B, L, and S denote the baryon number, the lepton 
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We consider the regeneration of gravitinos that starts after completion of reheating. The 
calculation of the production rate of these thermally produced gravitinos requires a con- 
sistent finite-temperature approach. A result that is independent of arbitrary cutoffs 
was derived for supersymmetric quantum chromo dynamics (QCD) in Ref. [2] . Following 
this approach, we provide for the first time the complete Standard Model gauge group 
result to leading order in the gauge couplings. For gravitino dark matter scenarios, 
this allows us to calculate the relic density of thermally produced gravitinos, ■ The 
comparison of this density with the observed dark matter density is crucial to decide 
whether gravitinos can explain dark matter. 

If the gravitino is not the LSP, it can decay at late times. The late decays of thermally 
produced gravitinos can then affect the abundances of light elements during big bang 
nucleosynthesis (BBN). Thus, the new results will also serve as central input parameters 
for deriving cosmological constraints for scenarios with unstable gravitinos. 

Thermal gravitino production becomes very efficient if the reheating temperature Tr of 
the Universe after inflation is high. Thus, gravitinos play an important role in models 
where thermal leptogenesis explains the matter-antimatter asymmetry of the Universe. 
At the dawn of the Large Hadron Collider era, we face the exciting possiblity to con- 
firm SUSY directly in experiments. In fact, we show in this thesis that a conceivable 
determination of the gravitino mass at future colliders will allow for a unique test of the 
viability of thermal leptogenesis in the laboratory. 

This thesis is organized as follows. Chapter 2 discusses the general supergravity La- 
grangian. We show how to obtain an effective low-energy theory. This allows us to derive 
the complete set of Feynman rules that is necessary for the subsequent calculations. In 
Chapter 3 we identify all processes for thermal gravitino production to leading order in 
the Standard Model gauge couplings. This yields the complete SU(3)c x SU(2)l x U(1)y 
result for the thermal gravitino production rate. Chapter 4 covers phenomenological 
implications of the new result for cosmology. In particular, for gravitino dark matter 
scenarios, we provide and derive an upper limit on the reheating temperature of 
the Universe. In Chapter 5 we propose the collider test of the viability of thermal lep- 
togenesis. Here, we take into account also the non-thermal gravitino production from 
decays of the next-to-lightest supersymmetric particle (NLSP) into the gravitino LSP. 
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Chapter 2 



Prom Supergravity to Gravitino 
Phenomenology 

The natural scale of supergravity is the Planck scale. For our purposes, however, we 
need to find a low-energy version which is reconcilable with phenomenology. On this 
quest we shall keep track of the gravitino interactions. This chapter therefore provides 
a trail from supergravity to gravitino phenomenology. 

After briefly sketching how the gravitino emerges in supersymmetric theories, we quote 
the general form of the supergravity Lagrangian^ in four spacetime dimensions and 
discuss its properties. We explain the conditions for spontaneous supersymmetry break- 
ing. For a concrete model of gravity-mediated supersymmetry breaking, we carry out 
the transition to the low-energy effective theory relevant for gravitino phenomenology. 
We relate the low-energy supergravity Lagrangian with the Minimal Supersymmetric 
Standard Model (MSSM) in the high-energy limit of unbroken electroweak symmetry. 
Having identified the couplings of the gravitino to the fields of the MSSM, we provide 
the complete set of Feynman rules necessary for the calculations in this thesis. 



We only consider unextended A'' = 1 supersymmetry. 



3 



2.1. Local Supersymmetry 



2.1 Local Supersymmetry 

Supersymmetry is a spacetime symmetry. It relates the bosonic and fermionic degrees 
of freedom in a supermultiplet of particles: 

Q I boson) ~ |fermion), Qjfermion) ~ jboson) , (2-1) 

The anti-commuting Weyl spinor Q, which is the generator of the SUSY transformation, 
extends the Poincare algebra to a graded Lie algebra. The generators Q fulfill the 
fundamental commutator [ • , • ] and anti-commutator { • , • } relations^ 

{Qa,Qp} = 2a>'^^P^, (2.2a) 
{Qa,Q(i} = {Qa,Qp} = 0, (2.2b) 

[P^,Qa] = [Pf.,Qp]=0. (2.2c) 

Let e and r] be two infinitesimal Weyl spinors which parameterize the supersymmetry 
transformations. We then can write (2.2a) in terms of a commutator: 

[eQ,fjQ]=2ea''fjPf, . (2.3) 

In a globally supersymmetric theory, the parameters e and i] are spacetime independent, 
corresponding to a rigid translation. Gauging supersymmetry corresponds to making 
the supersymmetry parameters local, i.e., dependent on x. Then, the commutator of 
two local supersymmetry transformations yields translations ~ e{x)a^r]{x) which differ 
from point to point. Invariance under these general coordinate transformations is exactly 
what we expect from a theory of gravity. Local supersymmetry is therefore referred to 
as supergravity (SUGRA)^. 

The gravitino is a central element of SUGRA because it is the gauge field of local super- 
symmetry transformations. Since supergravity is a non-renormalizable theory, operators 
of mass dimension five and higher occur. The gravitino thus is an extremely weakly 
interacting particle with couplings suppressed by inverse powers of the reduced Planck 
mass 

Mp = ^ = 2.4 X 10^^ GeV , (2.4) 
where Gn is Newton's gravitational constant [5]. 



^The conventions on spinor notation are found in Appendix A 

^Good reviews about supergravity are [3] and [4]. For further reading see references therein. 
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2.2. The Supergravity Lagrangian 



2.2 The Supergravity Lagrangian 

Our starting point of our analysis is the Lagrangian given in Appendix G in the book of 
Wess and Bagger [6]. We rewrite it in terms of four-component spinors and adopt it to 
our conventions given in Appendix A. This involves a change of signature in the Lorentz 
metric from (— , +, +, +), used by Wess and Bagger, to the signature more commonly 
used in high-energy physics (-I-, — , — , — ). With our conventions, we can carry out this 
task in a safe way. We explicitly restore units of Mp so that the supergravity Lagrangian 
in the four-component formalism then reads: 




+ - V2gdiDaX''x'L + \V2g [(Re/)"!] d^h^D^'xi 



8Mp 



mn 



la 



AMI 



1 



+hDiDjW^x{ + \g''* D,*W*difabX'}iX'i + h.c. 




(2.5) 
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2.2. The Supergravity Lagrangian 



Let us briefly introduce the building blocks of this Lagrangian: 

Matter sector Matter fermions are described in terms of left-handed four-spinors 

Xl = (^^^"l^^/") , (2.6) 

where (Xa)w/B stands for the two-component Weyl spinor. Here and in the fol- 
lowing, the subscript W/B indicates the quantities used in the book by Wess and 
Bagger [6]. Since the matter sector of the supergravity Lagrangian is built up 
in terms of left-chiral superfields, only left-handed four-spinors Xl appear in the 
resulting Lagrangian in the four-component notation. The corresponding scalar 
superpartners are denoted as The index i runs over all chiral superfields. 

Gauge sector The gauge multiplet consists of gauge bosons and their superpart- 
ners, the gauginos 

= I I , (2.7) 

V ^(-^''")w/B/ ' 

which are Majorana fields. Both are in the adjoint representation of the gauge 
group, a,b, - ■ ■ = 1, . . . , dim^. The associated field strengths of the gauge fields 
are written as F^^. The auxiliary fields Da are generalizations of the D-terms 
in the vector supermultiplets of a globally supersymmetric theory. 

Gravity sector The graviton e^™" shows up implicitly as the determinant of the vielbein 
e = det e^'" and the curvature scalar R. Flat spacetime indices of the local Lorentz 
frame are denoted by m, n, . . . ; Einstein indices by ^, z^, . . . . 

The gravitino is a spin-3/2 field which is written in terms of the Majorana vector- 
spinor, 

, /'-«(^m")w/b\ 

V hV'/. )w/b/ 

Since it is the superpartner of the graviton, the gravitino is massless for exact local 
supersymmetry. Note, that we include factors of i in the definition of the gravitino 
as well as for the gauginos in (2.7). 

Gauge kinetic function The gauge kinetic function fab{4>) is a dimensionless analytic 
function in the scalars (p. In the superfield approach, it multiplies the kinetic term 
for the vector supermultiplet. In the component version (2.5) of the supergravity 
Lagrangian, it therefore shows up as a prefactor in the kinetic terms of the gauginos 
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2.2. The Supergravity Lagrangian 



and gauge bosons, 

'-RefabX^rd^X' , (2.9a) 
-^Re^bF;,F^''^^ (2.9b) 

respectively. Because of gauge invariance, it transforms under gauge transforma- 
tions as the symmetric product of adjoint representations of the gauge group Q. 

Note that derivatives of the dimensionless gauge kinetic function fabi4>) with re- 
spect to the scalars 

have negative mass dimension. In the course of spontaneous symmetry breaking, 
it is reasonable to consider models where difab = O(Mp^). Therefore, we have 
not explicitly written out terms in the Lagrangian (2.5) which are proportional to 
oc Mp^ difab since they are of order 0(Mp ^). 

Superpotential The superpotential is an analytic function in the chiral scalar fields. 
It has mass dimension three and its general from is restricted by gauge invariance, 

W = VFh(/i) + ^/ii,(/i)<^V^' + ^y.,fc(/i)0V^</>' + 0(Mp 1) . (2.11) 

Here, we distinguish between the observable part of the superpotential and the 
hidden superpotential Wii{h). The latter depends only on hidden scalar fields h, 
which have no or only very small couplings to ordinary matter and gauge fields 
of the observable sector. The simplest choice is a superpotential with separated 
hidden and observable sectors 

W = Wi,{h) + Wo{^). (2.12) 

Let us remark here that a hidden field dependence f^ij{h) which mixes observable 
and hidden sector can be crucial in models offering a dynamical solution to the 
fj, problem [7, 8]. 

Kahler structure The Kahler potential K{(f>, (j)*) is a real-valued function in the chiral 
scalar fields and has mass dimension two. It has the generic form: 

K = Ki,{h, h*) + aij{h, h*)(j)'<t)*' + [Z,j{h, /i*)^^^' + h.c] + 0{Mp^) , (2.13) 

where i^h denotes the hidden part of the Kahler potential that is independent of 
observable fields. 
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2.2. The Supergravity Lagrangian 



In fact, the theory is endowed with a rich Kahler structure. One can think of the 
scalars as coordinates of a Kahler manifold whose metric gij* can be expressed 
by the second derivatives of the Kahler potential K{(j), (f>*): 

with its inverse g^^ , i.e., g^^ gj*k = The Kahler connection is given by 

r& = /'-^. (2.15) 

For example, the Kahler metric enters in the chiral sector in terms of a prefactor 
for the kinetic terms of the chiral fermions 

imrxii^d^^xl ■ (2.16) 

The Kahler covariant derivatives of the superpotential read 

DiW = Wi + Mp^KiW , (2.17a) 
ViDjW = Wij + Mp^{KijW + KiDjW + KjD^W) 

- T%DkW + 0(Mp 3) , (2.17b) 

where Ki = diK, Wi = diW, Kij = didjK, and Wij = didjW. Beside being 
invariant under local supersymmetry- and gauge transformations, (2.5) is inert 
under Kahler transformations with arbitrary holomorphic functions F{<j)) 

K{(l), r ) ^ Jf(0, r ) + F{^) + F*{^*) , (2.18a) 
^ e-^(<^)/*-^V , (2.18b) 

provided that the spinor fields undergo F-dependend Weyl rotations 

xl-e^i-^A^p^^X^^, (2.19a) 

^^-praF/M^jsx'^ ^ (2.19b) 

^^^e-i^-^/^^p^^^^. (2.19c) 

It is easily seen that the metric gij* and the Kahler connection T^j remain un- 
changed under a Kahler transformation (2.18). More generally, the isometries of 
the Kahler manifold can be expressed by the real Killing potentials Da{(j), 4>*)- The 
corresponding (holomorphic) Killing vector fields are given by 

= -ig'^*dj.Da , (2.20a) 
X*J = ig'^*diDa . (2.20b) 
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2.2. The Supergravity Lagrangian 



In the Lagrangian (2.5) we used this to express the Kihing vector fields in terms 
of the Killing potentials. In principle, the Killing potentials Da can be found by 
solving the corresponding Killing equation. We quote their explicit form later 
in (2.35), where we restrict ourselves to flat Kahler manifolds. 

The physical fields A'^ and A" in the vector multiplet are written with upper gauge 
indices a,b, . . . while the Killing potentials Da have lower indices. To bring them in 
canonical form, they are lowered or raised with Re fab or [(Re/)"^]'*'' , respectively. 

Furthermore, the supergravity Lagrangian above is given for a simple gauge group only. 
An extension to a non-simple gauge group Q = Ha like SU(3)c x SU(2)l x U(1)y 
requires to introduce an additional index. We will do so in Sec. 2.6. 

The covariant derivatives are defined as^ 

= d^cj)' + igg'^'dj.DaA''^, (2.21a) 

+ igd, [g''*dj,Da) Alxi + 0{Mp'), (2.21b) 
V^X" = d^X' - gr'^A'^X' + 0(Mp 2), (2.21c) 
V^iP, = df,ij, + 0{M-^). (2.21d) 

The Lagrangian (2.5) is invariant under local supersymmetry transformations, 
parametrized by means of an anticommuting Major ana spinor 



C(x) = ( (^(^^"^W/B 



(2.22) 



of mass dimension —1/2, namely. 



Jc^r = ^^^(^^'"^/^^ - ^^"^Mi) (2-23a) 

d^cP' = V2Cxl (2.23b) 

- V2e^^/'^^^pg'^'Dj,W*CL + 0{Mp^) (2.23c) 

5cA«=C7^A?i-C7MA2 (2.23d) 

•^cA^ = If^Ai^^YKl - ig [{Refr'Y' DaCi + 0{M-^) (2.23e) 

= 2MpP^C + ]^^e^^2A/2(^^^^^ ^ W*^^Cl) + 0{Mp^) (2.23f) 



*We drop the contributions from the spin-connection a;^™" in the covariant derivatives of the fermion 
fields since we wiU only consider flat spactime in this thesis for which Ui_J^^ 0. 
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2.3. SUPERSYMMETRY BREAKING 



Let us stress that the Lagrangian (2.5) depends on two arbitrary functions of the chiral 
scalars, namely, the Kahler function 

G(^,^.)= (2,24) 

and the gauge kinetic function /at (</>)• This is a consequence of the Kahler- Weyl invari- 
ance and can be seen by virtue of (2.18) with F{(j)) = \iiW{(j)). Within the framework of 
supergravity, there is no mechanism which tells us one particular form of G and fab from 
the other. On the other hand, phenomenological considerations will suggest a certain 
minimal choice as described below. 



2.3 Supersymmetry Breaking 

Particles within the same supermultiplet are degenerate in mass, because 

[P',Q^] = [P^Q^]=0. (2.25) 

Since we do not yet have experimental evidence for supersymmetry, we know that it has 
to be a broken symmetry if realized in nature. 

The condition for the spontaneous breakdown of a symmetry is a ground state that 
does not respect this symmetry. In the case of local supersymmetry, we find from the 
supersymmetry transformations (2.23) that the only Lorentz invariant ways to achieve 
this are^: 

{^CXl) Of {F' )C / F-term breaking (2.26a) 

and/or 

((5^ A") oc {Da)C + D-term breaking. (2.26b) 

Here, ( • ) denotes the vacuum expectation value (VEV) and W" is given by 

pi ^ e^''^^^g'^''Dj.W*. (2.27) 

This is a generalization of the auxiliary fields pf^'^^ = — c?i*M^* which are part of the 
chiral supermultiplets in a globally supersymmetric theory. 

Many models of spontaneous symmetry breaking have been proposed. Common to 
these models is that supersymmetry breaking occurs in the hidden sector of particles. 

^We do not consider the possibility of fermion-antifermion condensates which require the presence a 
strong gauge coupling force. 
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2.4. Gravity-Mediated Supersymmetry Breaking 



Super symmetry breaking is then communicated to the observable sector either at tree 
level or radiatively. In the course of spontaneous SUSY breaking, the gravitino acquires 
a mass mg. Depending on the SUSY breaking scheme, tuq can range from the eV scale 
up to scales beyond the TeV region [9]. 

2.4 Gravity-Mediated Supersymmetry Breaking 

An appealing scenario of supersymmetry breaking is gravity-mediation. We show how 
phenomenological considerations can lead to certain minimal choices of the gauge kinetic 
function and the Kiihler potential in such a scenario. Furthermore, we illustrate the 
super-Higgs mechanism and show how soft super symmetry-breaking terms arise in the 
low-energy limit of supergravity. Finally, we also comment on the way Yukawa couplings 
emerge. 

Supersymmetry is broken by VEVs (F™ ) of some scalar fields h which are gauge singlets. 
The hidden fields h appear in non-renormalizable terms of mass dimension higher than 
four, which are suppressed by powers of Mp. We expect soft terms of order 



to arise. This is intuitive since mgoft has to vanish in the limiting case of unbroken super- 
symmetry (F) 0. In order to obtain msoft in the electroweak range, supersymmetry 
breaking has to occur on r ather high scales, ^/{F) ~ 0(10^° GeV). 




In the following we will sketch such a scenario and derive a low-energy limit of the 
supergravity Lagrangian (2.5). The phenomenological imperative for any concrete 
supersymmetry-breaking scenario is that it delivers a supersymmetric theory equipped 
with terms that softly break it. The soft terms obtained from our chosen gravity media- 
tion scenario will have specific sizes which might differ significantly from other breaking 
scenarios. Nevertheless, their structure will be the same. In the end, we therefore can 
treat the soft parameters as free parameters (to be determined experimentally) and for- 
get about the details of the underlying breaking scenario. What we will have gained is 
a consistent set of Feynman rules for the MSSM in the high-energy limit together with 
a set of gravitino rules which share the same convention. 

Let us start with the gauge kinetic function. The trivial choice fab = ^ab would render 
the kinetic term of the gauginos (2.9a) and the kinetic term of the gauge bosons (2.9b) 
renormalizable. This form, however, is phenomenologically disfavored since we cannot 
produce gaugino masses at the tree level. ^ The candidate for the gaugino mass term in 
®A tree-level gaugino mass can also arise from D-term breaking which we do not consider in this 



"^soft 





Mp 
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2.4. Gravity-Mediated Supersymmetry Breaking 



(2.5) is 

le'^/^^'h''* D,,W*difaiA'kXl + h.c. , (2.29) 

which vanishes for fab = Sab- We can obtain a diagonal gaugino mass matrix with the 
generic ansatz: 

fab = 6abf{h) . (2.30) 

After some hidden field h has aquired a VEV {h) at some high scale, we obtain canonically 
normalized gauge kinetic terms by a rescaling 



A*^ = v'(Re/)A" (2.31a) 
A1 = vW)^;i- (2.31b) 



Rewriting (2.5) in terms of A" and A'^, we see immediately that 1/-^/ (Re/) appears in the 
covariant derivatives as prefactor of the gauge coupling g. This is remarkable because it 
suggests that Re/ plays the role of a coupling constant. The rescaling 

then completely hides any dependence on He fab in (2.5). If we further assume that (h) is 
real, all terms proportional to luifab vanish in (2.5). For our purposes, this assumption 
is safe since only terms of 0{Mp^) with a gravitino in the vertex will be considered; and 
an operator cx Mp^lmfabipfi is absent in the Lagrangian (2.5). 

The term (2.29) contains also the superpotential W and the Kahler potential K. The 
general form of the Kahler potential is specified in (2.13). There, a diagonal choice 
aij{h, h*) = a{h, h*)6ij is phenomenologically favored. A non-diagonal form of Oij leads 
to off-diagonal terms in the the scalar mass matrix in the low-energy limit. This can 
induce flavor-changing neutral currents for which strong bounds exist [9]. 

In the case of the MSSM emerging as a low-energy effective theory, we can have a con- 
tribution like [Z{h, h*)HiH2 + h.c] in the Kahler potential where Hi and H2 denote the 
the two SU{2)l Higgs doublets. The hidden field dependence of Z is closely connected 
to the fj, problem [7, 8], but is beyond the scope of this thesis. A particularly appealing 
choice of the Kahler potential is therefore given by 

K = Kh{h, h*) + Y^ + [Z{h, h*)HiH2 + h.c] (2.33) 

i 

with aij = 6ij [see Eqn. (2.13)], corresponding to a flat Kahler manifold in the observable 
sector. With this choice, the Kahler metric (2.14) becomes trivial, 

9ij* = ^ij* ; (2.34) 

work. 
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2.4. Gravity-Mediated Supersymmetry Breaking 



and leads to canonical kinetic terms in (2.16). Moreover, the connection coeffi- 
cients (2.15) vanish: F^^ = 0. The Killing potentials Da then coincide with the D-terms 
of a globally supersymmetric Yang-Mills theory: 

Da = n,ij^- (2.35) 

The generators of the Lie algebra of the gauge group are denoted by T^. They are chosen 
to be hermitian and do obey the commutator relation 

[Ta,n] = ir'^T, (2.36) 

with the structure constants f"-'^'^. 
Super-Higgs mechanism 

In supergravity, an analogon to the Higgs mechanism of electroweak-symmetry breaking 
exists. When supergravity is spontaneously broken, the corresponding massless Gold- 
stone fermion, or goldstino, is absorbed by the gravitino which aquires thereby its ±1/2 
helicitiy components. 

Mass terms for the gravitino, i.e., terms in (2.5) involving the gravitino field, which are 
quadratic in the fermionic fields and do not have derivative couplings are 

- ^e^/'''''^W*4^j,^[r,Y]i^Lu - ^e^/2^^^p V2Al^V^^7'^xi + h.c. . (2.37) 

Since hidden sector fields do not share any gauge interactions, there are no D-terms 
contributing to supersymmetry breaking via (Da) / 0. Therefore, we do not consider 

^ -gDaX1i-/''^l^f, + h.c. (2.38) 



2Mp- 

as a potential mass term. 

Because of their hidden field dependence, the gauge kinetic function, the superpotential, 
and Kahler potential will get vacuum expectation values (fab), (I^h), and (Wh), respec- 
tively. From the first term in (2.37) we can read off the gravitino mass after spontaneous 
symmetry breaking: 

ma = -^^^""'^^'''^^ (W^) . (2.39) 

The second term in (2.37) mixes the gravitino with chiral fermions Xl- Let us define 
the spinor^ 

VL = DiW^L (2.40) 



^In case of D-term breaking, the definition (2.40) of the goldstino has to be extended in order to 
remove the mixing in (2.38). 
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2.4. Gravity-Mediated Supersymmetry Breaking 



and apply the supersymmetry transformation (2.23). We see that t] changes by a shift 
S^rjL = -V2e''/^'h''* iD^W) {D,,W*) + ■■■ = -3V2m^C + ■■■ , (2-41) 

where in the last equality we have assumed a vanishing cosmological constant as ex- 
plained below. This implies that 

r? = AWxi + A*^^*(xi)' (2.42) 

is the Goldstone fermion. Indeed, with the choice ^ = r/\/2/(6mg), we can choose a 
unitary gauge where rj transforms to zero. This removes the mixing term in (2.37). 

Soft terms from the scalar potential 

The scalar potential in the supergravity Lagrangian (2.5) reads 



iPiW) (Z),-.VF*)-3 



Ml 



= g.,.rF^* - 3e^/^i ^ + '-g^D^D^ , (2.43) 

where the rescaled coupling (2.32) is used. We have already seen that successful symme- 
try breaking is achieved if some of the auxiliary fields (2.27) get a non-zero vacuum 
expectation value 

(F') = (e^l^^^^g^^*Dj,W*^ / . (2.44) 

While the scalar potential 

^giob = FF* + l/2D^ (2.45) 

is positive semi-definite in global supersymmetry, the minus sign of the second term in 
(2.43) offers the appealing possibility to break supersymmetry with vanishing cosmolog- 
ical constant: 

{V) = (g^^* {D,W) {D,.W*) - 3^^ = . (2.46) 

We therefore can work in the limit of flat Minkowski spacetime and neglect the inter- 
actions of the graviton: — > and e — > 1. Moreover, we drop the distinction between 
flat and curved spacetime indices. Because of notational habit, we denote from now on 
Minkowski spacetime indices with /x, = 0, . . . , 3. 

For the sake of simplicity, we further assume that we do not have any mixing between 
the observable and hidden sectors in the superpotential and the Kahler potential. For 
the superpotential, this corresponds to the choice (2.12). For the Kahler potential, we 
write 

K = K^{Kh*) + Y,<t^'^*' . (2.47) 
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2.4. Gravity-Mediated Supersymmetry Breaking 



It is instructive to get a feeling for the scales involved: 

{h) - 0(Mp), (P^h) ~ 0(M2), (K^) ~ 0(m2), (a^W'h 



0{Mp). (2.48) 



The second and third relations come from the requirement that we want to obtain a 
gravitino mass (2.39) in the electroweak range < few TeV. This requirement is 
crucial because rriQ will govern the size of the soft terms. The fourth relation is a 
consequence of the second. 

In the following we will use indices m,n . . . for fields h which we allow to aquire a VEV 
and indices for observable fields (j). Using the definitions (2.17) for the Kahler 

covariant derivatives, the scalar potential (2.43) without the D-term becomes 



Vf 



E 



+ 



m 



+ 



1 dK^ 



{Wo + w^) 



Ml 



(2.49) 



We expand this to 
Vf =e 



E 



+ 



M4 



+ 



^ dW* 



Ml 



+ E 



Ml dh*-^ dh"^ 



Ml a/i*™ 



' Ml 



' Ml 



+ 0{Ml) + 0{Mp^] 



(2.50) 



To get from (2.49) to (2.50), we have expanded the exponential of the observable sector 
in powers of Mp ^. 

Now we can the perform flat limit of supergravity where one sends the Planck mass Mp to 
infinity but holds the gravitino mass (2.39) fixed. Terms of 0{Ml) in (2.50) contain 
relative signs and have to be fine-tuned to zero for vanishing cosmological constant 
{V) = 0. With a rescaling of the observable superpotential. 



Wo = Wo 



(2.51) 



this yields 
Vf = 



E 



dWo 



+E 



Ml 



dh*" 




(2.52) 
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2.4. Gravity-Mediated Supersymmetry Breaking 



In the first term of (2.52) we recover tfie scalar potential of global supersymmetry (2.45), 
where Ff"^ = -di^W*. (We do not carry along the D-tevms which are already shown 
to coincide with the globally supersymmetric case in (2.35).) From (2.27), we find 

Let us identify the remaining parts in (2.52) with the commonly used notation^: 

mo = mg , (2.54a) 

In Sec. 2.6 we will write down the superpotential of the MSSM. It is then not hard to 
show that A corresponds to the trilinear scalar coupling and that i? is a bilinear mass 
parameter; mo is a universal scalar mass. The values of the above soft-supersymmetry 
breaking parameters are in the electroweak range. Since the gravitino mass governs 
the size of the soft terms, the choices made in (2.48) turn out to be phenomenologically 
favorable. 

Gaugino mass term 

To complete the set of soft breaking parameters (2.54), we discuss now the gaugino mass 
term (2.29). Following the logic above, it is straightforward to perform the flat limit. 
One then finds in terms of the rescaled gaugino fields (2.31a): 

1 



^MabXIiXl + h.c, (2.55) 



with the mass matrix [cf. (2.30)]: 



From (2.53) together with (2.48), one finds that (F™) ~ 0{Mp) so that the gaugino 
masses (2.56) are in the electroweak range for 

(Refih)) = 0{l) and ^M^^^ ^ ©(Mp^). (2.57) 

For example, a straightforward choice is f{h) = h/M-p. 



*The relation B — A — is just a consequence of the simple breaking scenario which we have 
chosen. In general, one treats A and B as independent parameters. 
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2.5. The Free Gravitino Field 



The second relation in (2.57) justifies our dropping of terms oc Alp^difab in tlie La- 
grangian (2.5), because tliese terms are of C'(Mp^). 

The flat limit, i.e., Mp — > oo and fixed, amounts to integrating out fields in the 
Lagrangian which are suppressed by powers of Mp. We therefore expect that the re- 
sulting Lagrangian is valid at a high energy scale. Thus, the soft parameters (2.54) and 
(2.56) should be interpreted as boundary conditions for the renormalization group (RG) 
equations [10]. 

Yukawa interactions and chiral fermion mass terms 

As far as spontaneous supersymmetry breaking is concerned, there is another term in 
(2.5) to address, namely, 

AeK/2M^.V,D,W^xi + h.c. . (2.58) 

The definition for the Kahler covariant derivatives is given in (2.17). Recalling the 
dimensional considerations (2.48), one finds that the only contribution for the observable 
sector fields that survives in the flat limit is ^ 

1 d'^Wo i i . , . 

Here we made use of the rescaling of the observable superpotential (2.51). This rescaling 
can be absorbed into the parameters of the superpotential (2.11), 

We see that the bilinear part of the superpotential (2.11) produces mass terms for chi- 
ral fermion flelds with mass parameter fiij while Yukawa couplings yijk arise from the 
trilinear part of the observable superpotential. 



2.5 The Free Gravitino Field 



In the preceding section we have seen that the gravitino acquires a mass through 
the super-Higgs mechanism. From (2.5) together with (2.39), we find the Lagrangian for 
the free gravitino field 

jrtee ^ e^^P'^i^i^^-f^dp^L. " -m^ (?i^^[7^ tIV'l. + h.c.) (2.61) 
^In fact, a Weyl rotation (2.19a) of the fermion fields with F — — -^^y^ — Afp In ), is necessary. 
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2.5. The Free Gravitino Field 



which can be rewritten as 



^free ^ J_^,upa^^^^^^Q^^^ _ ^ ^^^^^ [y^ , ^-]^^ + tot.div. . (2.62) 



Variation of (2.62) yields the Rarita-Schwinger equation [11] 



l^ludpi^a - -:rnQ[j'',Y]i^, = . (2.63) 



Since the gravitino satisfies the constraints 

T^V'm = , (2.64a) 
a^V'M = , (2.64b) 

the Rarita-Schwinger equation (2.63) can be shown to reduce to the Dirac equation for 
each vector component /x of the gravitino, 

{i$ - mg) = . (2.65) 

The polarization tensor for a gravitino with four-momentum P is given by [2] 

s 

(2.66) 

where the sum is performed over the four gravitino helicities s = ±3/2, ±1/2. The 
polarization tensor obeys 

^^n^u{P) = , (2.67a) 

P^U^uiP) = , (2.67b) 

(i^-m)n^^(P) = . (2.67c) 

For energies much higher than the gravitino mass, it can be shown that the II^i, splits 
into two parts [2], 



3 m?- 
G 



n^,(P) ~ -fg^, + -f^^ . (2.68) 



The first term represents the sum over the helicity ±3/2 states of the gravitino whereas 
the second part corresponds to the sum over the ±1/2 helicities of the goldstino (2.40). 
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2.6. The MSSM in the High-Energy Limit 



2.6 The MSSM in the High-Energy Limit 

In the calculation of the gravitino production rate we will assume a primordial 
plasma with the particle content of the Minimal Supersymmetric Standard Model 
(MSSM). Since we will consider thermal production far above the electroweak symmetry- 
breaking scale 0(100 GeV), we can work with an unbroken electroweak symmetry group 
SU(2)l X U(1)y. 

We now collect our previous results and extend the above considerations to the standard 
model gauge group, i.e., we introduce an additional index a that keeps track of the 
different factors 

3 

a = n = U(1)y X SU(2)l X SU(3)c (2.69) 

SO that henceforth we use the assignment 

a = I for U(1)y , a = 2 for SU(2)l , a = 3 for SU(3)c . (2.70) 

Accordingly, the gauge couplings are given by 

9i=g' , 92 = 9 , 93 = 9s , (2.71) 

with the U(1)y hypercharge coupling g' , and the weak and strong coupling constants g 
and gs, respectively. 

In terms of the rescaled quantities (2.31) and (2.32), we find from (2.5) for the SU(3)c x 
SU(2)l X U(1)y gauge interactions 

3 

■^gauge = -^gauge (2.72) 
a=l 

with 



4aige = +-Di^U"D^'''>'r - ^9l{rT^l^'y 



- V25aA(") >«ri") xi - V2g^tLTi% A(°) ^ . (2.73) 



Note that we have dropped the hats which were introduced originally to indicate the 
rescalings, i.e., we set 

g^g, A^A, A^A. (2.74) 



19 



2.6. The MSSM in the High-Energy Limit 





Table 2.1: Gauge 


fields of tlie MSSM 








1-1 a 1 1 cci n ot; A v*^/ ^ 




B-boson, bino 






(l,l)o 


W-bosons, winos 






(l,3)o 


gluon, gluino 


A (3) a _ 


A(3) a —ga 


(8,l)o 



With the D-term (2.35) and the vanishing Kahler connection for trivial Kahler manifolds, 
the covariant derivatives (2.21) in the Lagrangian (2.73) become 



The field strength tensor F^^^ " reads 



o) c 



(2.75a) 
(2.75b) 
(2.75c) 



(2.76) 



The corresponding gauge fields and their superpartners are listed in Table 2.6. Matter 
fields with gauge couplings are in the fundamental (anti-fundamental) representation 
of the corresponding gauge group, namely, 2 (2 = 2) for SU(2)l and 3 (3) for SU(3)c. 
The SU(2)l doublet structure for the MSSM matter fields is reviewed in Table 2.6. The 
strongly interacting particles gather in color triplets 3, i.e., squarks and quarks indicated 
in Table 2.6 carry an additional color index. Gauge singlets are denoted by 1 and (•, •)o 
for the strong/weak and U(1)y interactions, respectively. Note that the normalization 
for the hypercharges is such that the electric charge Q is given hj Q = + Y/2 where 
T3 denotes the weak isospin eigenvalue = ±1/2 for upper/lower entries in the doublets 
of Table 2.6, respectively; T3 = for SU(2)l singlets (•,l)y. The family index / in 
Table 2.6 refers to one out of three generations of, 



leptons < 



-cl 



( '^1 I^R ' 



sleptons < 



-R 
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2.6. The MSSM in the High-Energy Limit 



Table 2.2: Matter fields of tlie MSSM 



Name Bosons 0* Fermions xl (SU(3)c, SU(2)l) 



Y 



I - [ "^L \ tI 



Sleptons, leptons = ( _^ J = _^ J ( 1 , 2 )_i 

1=1,2,3 



E*^ = ei*' E-' = el'' (1,1) 



+2 



^^^^^ 

(x 3 colors) = [/c/ = ^c/ ( 3 , 1 )_ 

D*^ = d*j( D^^ = d%^ (3,1)_ 



Higgs, higgsinos Hd={t\ = I 1 ( 1 , 2 ) _i 



( Ui , CL , tL ) 

/ - — ^5): - — if: ^* "\ 

{dL ,SL ,bL ) 

for I = (1,2,3), respectively. 

All matter fields are written in terms of left-handed four-spinors since they stem from left- 
chiral supermultiplets. For example, a right handed tau lepton with hypercharge —2 
is written in terms of its charge conjugate r^'^ with hypercharge +2. 

The generators in (2.73) for the standard model gauge group read 

Ta% = \yAM for U(1)y , (2.77a) 

Ti%=l^a,^J forSU(2)L, (2.77b) 

Ti%=lK^j forSU(3)c, (2.77c) 

where Yi is the hypercharge of the corresponding particle or x^l (^^^ Table 2.6). 
The Pauli sigma matrices aa (a = 1,2,3) are given in (A. 7) in the Appendix and 
Aa (a = 1, . . . , 8) denote the Gell-Mann matrices. The chosen basis (2.77) for the gener- 
ators implies that the structure constants /^"^ "'"^ for the non-abelian groups are totally 



quarks < 



u 



cl_ 



di- 



[di ,SL,bL) 



squarks < 



di- 
d*R'- 
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2.7. Gravitino Interactions 



antisymmetric. Since U(1)y is abelian, the commutator relation (2.36) is trivial, i.e., 

fWabc-Q forU(l)Y. (2.78) 

The superpotential of the MSSM is given by 

Wo = H^MSSM = U* YuQ -Hu-D* yaQ -Ha-E* y,L • + fiH^ • . (2.79) 

The doublet structure is tied together as Q ■ = e'^^QiHuj, with e*-' given in (A. 4a). 
Furthermore, U* yuQ is meant to be a matrix multiplication in family space, U* yuQ = 

The couplings yu,d,e and the bilinear parameter fi in (2.79) are understood to be the 
rescaled quantities on the left-hand sides in (2.60). We see from (2.59) that yu,d,e yield 
Yukawa coupings ~ y(t>Xi,XL and that /i gives mass to the higgsinos. Neglecting all 
Yukawa couplings except the one for the top quark, 

/O o\ 

ye ^ , yd ^ , (y^) ~ 




VO ytj 



(2.80) 



the superpotential (2.79) reduces to 

T^MSSM ^ ytt*RtLHl - ytt*nhLHt + ^ {HtHl, - H^H^^) . (2.81) 

In this thesis, however, we will not consider contributions to the thermal gravitino pro- 
duction rate coming from top- Yukawa interactions; we leave such an analysis for future 
work. 



2.7 Gravitino Interactions 

Let us now collect the terms in the supergravity Lagrangian (2.5) which describe the 
gravitino interactions with ordinary matter fields. In the calculation of the gravitino pro- 
duction rate, we will only consider external gravitinos which are subject to the constraint 
(2.64a). Therefore, the relevant interaction Lagrangian extracted from (2.5) is 



" ^^m[^''^'']^^^^"^''^p"^" • (2.82) 

Here we have used the Kahler metric (2.34) and the definition for the rescaled 
fields (2.31), but dropping the hats as in (2.73). Note that each operator in (2.82) 
is suppressed by Mp^. 
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2.8. Effective Theory for light Gravitinos 



2.8 Effective Theory for light Gravitinos 

In the discussion of the super-Higgs mechanism, we have seen that the goldstino degrees 
of freedom become the hehcity ±1/2 components of the gravitino. Indeed, these longi- 
tudinal helicity ±1/2 components become dominant if the energies involved are much 
higher than the gravitino mass mg. The dynamics of the Goldstone fermion is given by 
the derivative coupling of the goldstino to the supercurrent. The correct effective La- 
grangian for light gravitinos with goldstino-matter couplings in non- derivative form has 
been found in [12]. For a single external goldstino the non-derivative form is equivalent 
to the derivative form to all orders in perturbation theory and reads 



9a Ma 



rri°UV75A(-)V (2.83) 



where the Major ana goldstino field is denoted by ijj}^ Note that all vertices are pro- 
portional to supersymmetry-breaking mass terms. The coupling in the first term is 
proportional to the squared masses m^^ and m^; of the corresponding matter fields 0* 
and x\- The couplings in the remaining terms are linear in the gaugino masses M„. 
Therefore, at high energies and temperatures T, contributions involving the goldstino- 
fermion-scalar vertex are suppressed relative to the gaugino contributions due to the 
higher mass dimension of the coupling. 



2.9 Feynman Rules 

We are now in a position to provide all Feynman rules necessary for the calculation 
of the thermal gravitino production rate. The gauginos A", the gravitino V/xi and the 
goldstino "0 are Majorana fermions. Since these fields are self-conjugate, they yield 
different Wick contractions from those of Dirac fields. We therefore use the method 
proposed in [13] and introduce a continuous fermion flow, i.e., an arbitrary orientation 
of each fermion line. Proceeding against the fermion flow then allows one to form chains 
of Dirac matrices such that the relative sign of interfering diagrams can be obtained 
in the same manner as one does for Dirac fermions. We therefore have two analytical 
expressions for each vertex corresponding to the two possible orientations of the fermion 
flow. 

^"The definitions for the gaugino fields in [12] already contain the factors of i of our convention (2.7). 
For the goldstino tp, we have included them in the transition to the four-component formalism. 
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2.9. Feynman Rules 



For one direction of the fermion flow, the gravitino Feynman rules derived from (2.82) 
are given in Fig. 2.0. The gravitino is represented as a double solid line, scalars are 
given by dashed lines and chiral fermions x\ ^-I's given by solid lines. Gauge bosons A^°''> " 
are shown as wiggled lines and the corresponding gauginos are depicted as wiggled lines 
with additional straight solid lines. All momenta are understood to flow into the vertex. 



6, p,P 




4M, 




Figure 2.1: Feynman rules for the gravitino from (2.82) for one direction of the 
fermion flow. The wiggled lines represent gauge fields ^'■"^ while the corresponding 
superpartners A*^"^" are depicted by wiggled lines with additional straight solid lines. 
Scalars chiral fermions and the gravitino are represented respectively by dashed, 
solid, and double-solid lines. All momenta flow into the vertex. 



Figure 2.1 shows the relevant Feynman rules for the gauge interactions of the gauge group 
Qa derived from (2.73) for one direction of the fermion flow. For U(1)y, or a = 1, there 
is no self-coupling of the gauge bosons A^^'^ = B^. Therefore, the triple gauge boson 
vertex and its supersymmetric counterpart in the second line of Fig. 2.1 are absent. 

For a light gravitino, its interactions are dominated by goldstino dynamics. The cor- 
responding Feynman rules derived from the effective Lagrangian (2.83) are shown in 
Fig. 2.2. The double solid lines now represents the Majorana field ip. 

The complete set of the relevant Feynman rules including vertices with both directions 
of the fermion flow are given in Appendix B. 
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2.9. Feynman Rules 



K 




J 

a, v, P 

+ gf"'{K" -Q-)] 



V2gaT^lPL 




iV2gaT^liPR 



Figure 2.2: Relevant Feynman rules for the supersymmetric gauge interactions derived 
from (2.73). As in Fig. 2.0, it is understood that the gauge bosons and gauginos are 
associated with the corresponding gauge group (a). All momenta are ingoing. 



b, p,P 




m% — mfj 



Sab [AtI 



c,p 



9aMa (a) 



9 9 
ml; — mZj 




Figure 2.3: Feynman rules for the effective theory of light gravitinos from (2.83) for 
one direction of the fermion flow. The double solid line now denotes the Majorana 
spinor ip. Again, all momenta are understood to be ingoing. 
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Chapter 3 



Thermal Gravitino Production 



Thermal gravitino production in a consistent thermal field theory approach has been 
worked out for super symmetric quantum chromodynamics in Ref. [2]. Taking into ac- 
count also the electroweak processes, we extend the calculation to the full standard 
model gauge group SU(3)c x SU(2)l x U(1)y to leading order in the gauge couplings. 
We also correct an error in the SU(3)c result of Ref. [2]. 



3.1 The Braaten— Yuan Prescription 

In the previous chapter we have seen that gravitino interactions are suppressed by inverse 
powers of Mp. Thus, the dominant contributions to gravitino production and annihila- 
tion processes to leading order in the gauge couplings are inelastic 2 — > 2 reactions with 
one external gravitino. For one of these 2 — > 2 scattering processes, the net gravitino 
production rate at finite temperature reads [14] 



dTg _ 1 f d^p 



d^p 2(27r)3E J An 



i (27r)32^i 



{27r)^5\Pi +P2-P3-P) 



X [ME,)h{E2)[l ± f3m][l - fam\M{l + 2^3 + G)\' 

-[1 ± ± f2{E2)mE3)fQ{E)\M{3 + G ^ 1 + 2)|2} , (3.1) 

where E and P are the energy and four-momentum of the gravitino, respectively. The 
corresponding squared matrix element |Mp is weighted with the phase space distribution 
functions fi{Ei) of the particles involved in the scattering; it applies for final-state 
bosons/fermions and corresponds to Bose enhancement/Pauli blocking, respectively. In 
expression (3.1), an average dQp/iir over the directions of the gravitino momentum is 
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3.1. The Braaten-Yuan Prescription 



taken. The squared matrix elements \M\^ are assumed to be summed over initial and 
final polarizations and to be weighted with the appropriate multiplicities . 

At high temperatures, all particles except the gravitino are in thermal equilibrium so 
that /i, /2, and /a are given by the equilibrium distributions, 

1 



fB{Ei) 



1 



1 



where T denotes the temperature of the thermal bath. 



for bosons, 
for fermions. 



(3.2a) 
(3.2b) 



Assuming that inflation governed the earliest moments of the Universe, any initial pop- 
ulation of gravitinos must be diluted away by the exponential expansion during the 
slow-roll phase. We consider the thermal production (or regeneration) of gravitinos that 
starts after completion of reheating at the temperature Tr. Accordingly, the gravitino 
phase space density is much smaller than the equilibrium distribution /p. We thus 
can set (l — /g) ~ 1 and neglect gravitino disappearance processes 3 + G ^ 1 + 2. The 
production rate (3.1) then becomes 

3 



1 



d^p 2{2ttYE J Att 



dUr 



n 



d^Pi 



{2n)^5\Pi + P2-P3-P) 



h/B{Ei)h/B{E2)[l ± h/B{Ez)]\M{l + 2^3 + G)f 



(3.3) 



Note that a naive application of perturbation theory can lead to logarithmically singular 
contributions to the production rate. If a massless gauge boson with three-momentum 
|k| = A; is exchanged in the t- or n-channel, the corresponding squared matrix element 
is divergent for A; — > 0. A rigorous method to deal with such situations is the Braaten- 
Yuan prescription [14] which requires the weak coupling limit, g <^ 1: One introduces 
an intermediate momentum scale k* such that gT <^ k* <^ T. This scale seperates soft 
gauge bosons with momentum transfers of order gT from hard ones with momentum 
transfers of order T. The gravitino production rate is then given by the sum 



dT^ 



d^p d^p 



+ 



hard 



d^p 



(3.4) 



soft 



The hard part is obtained conveniently by computing the squared matrix elements 
|M((l + 2 — > 3 + G)p in standard zero-temperature perturbation theory. Whenever 
a massless gauge boson is exchanged in the t- or li-channel, k* will be introduced as an 
infrared momentum cutoff. The resulting production rate will then be of the form 

k* 



d^G 



d^p 



d^G 



hard 



d'^p 



^hard + In 



k*<k 



(3.5) 



27 



3.2. Hard Contribution 



In the region of soft momentum transfer, k < k*, one employs the finite-temperature 
version of the optical theorem. At T = 0, the imaginary part of the self-energy is 
related to the decay width of an unstable particle. In a thermal plasma, however, even 
stable particles can disappear and be produced in inelastic scattering off particles in the 
thermal background. Accordingly, the discontinuity in the self-energy gives the rate at 
which a non-equilibrium distribution of the corresponding particle approaches thermal 
equilibrium [15]. The soft part of the thermal production (or regeneration) rate of 
gravitinos can thus be expressed in terms of the thermal gravitino self-energy 



^ 1 fF{E) 

soft 



^ImSg(S + ie,p)|fc<fc. . (3.6) 



In the naive consideration of the production rate, the singular behavior stems from long- 
range forces. The physical cutoff is provided by the screening of the interactions due 
to the cooperative motion of particles in the thermal bath [14]. Instead of using bare 
gauge boson propagators, hard thermal loop (HTL) resummed gauge boson propagators 
have to be used in the region of soft momentum transfer, k < k* [16]. For gravitino 
energies E > T, there is no need to consider HTL-resummed vertices. The soft part of 
the thermal gravitino production rate (3.6) will then be of the form 



A,,i, + Bln(—] (3.7) 

\mth. 



soft 

with mth denoting the thermal mass of the corresponding gauge boson. 

Indeed, in the sum (3.4) the logarithmic dependence on k* cancels out and yields a 
finite result. The Braaten-Yuan prescription [14] together with the HTL-resummation 
technique [16] allows us to calculate the gravitino production rate in a gauge-invariant 
way. 



3.2 Hard Contribution 



From the Feynman rules for the gravitino and the supersymmetric gauge interactions 
presented respectively in Figs. 2.0 and 2.1, we find the leading processes for thermal 
gravitino production. Figure 3.2 shows these processes for one factor (2.69) of the 
standard model gauge group SU(3)c x SU(2)l x U(1)y- 

Note that the processes A, B, and F are not present for a = 1, i.e., for U(1)y, since 
there is no self-coupling of gauge bosons for an abelian gauge group; thereby also the 
supersymmetrized version of the relevant three-gauge boson vertex is absent. The matter 
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3.2. Hard Contribution 



A (BBF): ^(°)'^ + ^(°)* ^ A(")^ + G {$ for U(1)y) 






B (BFB): " + A(") ^ A^") ^ + G (crossing of A, $ for U(1)y 
C (BBF): 0* + ^ xi + G 




G 




G 




G 



D (BFB): ^(") + xl ^ (A'' + G (crossing of C) 
E (BFB): + xi ^ A(") + G (crossing of C) 
F (FFF): A(") " + A^") ^ ^ A^") + G (^ for U(1)y 



G A'")" 



G (FFF): xi + A(")'^^xi + G 





G A("'° 



A'"'' 




A{a)a 

H (BFB): + A(") 






G 




G 



I (FFF): xi + Xl^ ^ A(") + G (crossing of G) 
J (BBF): (p^ + <p*^ A(°) + G (crossing of H) 



Figure 3.1: The 2^2 scattering processes for gravitino production. Processes A, B, 
and F are not present for U(1)y since fW"-'"^ = 0. 
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Table 3.1: Squared matrix elements for the 2-^2 scatterings in terms of the Mandel- 
stam variables s and t; Ma denote the gaugino masses. Processes A, B and F are not 
present for U(1)y. Sums over initial and final spins have been performed. 

Label i Class Process i I-A^ip/ ^ ( 1 + 



3m4 



A 


BBF 




B 


BFB 


_)_ 


C 


BBF 




D 


BFB 




E 


BFB 


+ xi 


F 


FFF 


_|_ xi»)b 


G 


FFF 




H 


BFB 




I 


FFF 


Xl + Xl 


J 


BBF 





-+G 2.|ii::)p 



, {s'^+st+t^f I f(a)ahc\2 



(") 12 



2(. + 2t + 2$) \T^l^ 



fields and Xl coupling to the gauge bosons A^"^^ ° of the corresponding group Qa can 
be read from Table 2.6. 

The squared matrix elements for the processes shown in Fig. 3.2 are given in Table 3.2 
in terms of the Mandelstam variables 

s = {Pi + P2? , (3.8a) 
t = (Pi - P^f , (3.8b) 

where the four-momenta Pi, P2, and P3 are associated with the particles in the order 
in which they are written down in the column "Process i" of Table 3.2. The gaugino 
masses are written as Mq. 

The Dirac traces occurring in the evaluation of the squared matrix elements have been 
performed using the computer program FORM [17]. Polarizations sums over final and 
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initial states have been carried out (no averaging).^ In process A, there are two gauge 
bosons with four-momenta Pi and P2 in the initial state. The simple replacement of 
the corresponding polarization sum Ylpol ^a^V^ ~^ ~9^^5ab would amount to the inclu- 
sion of longitudinal polarizations. Since we work in the high-energy limit of unbroken 
electroweak symmetry, not only the gluons but also the electroweak gauge bosons are 
massless. Therefore, we have used the proper polarization sum 

5]e^(Pi/2)er(A/2) = 

pol. 

in which the unphysical modes are already subtracted [18]. The Feynman rules used in 
the evaluation of the Feynman diagrams are given in Appendix B. 

The squared amplitudes shown in Table 3.2 have to be weighted with appropriate mul- 
tiplicities: 

Processes A, B, F The sum over the gauge- group indices a, 6, c = 1, . . . ,d\mQa has 
to be performed: 

J^|/(")"'l' = iVa(iVi-l) for SU(iV„). (3.10) 

a,b,c 

For processes A and F, an additional factor of 1/2 occurs because of identical 
particles in the initial state. Note that the structure constants vanish for U(1)y- 

Processes C, D, E, G, H In the corresponding initial state, there can be 

2n2 = 28, (3.11a) 

distinct isospin-doublets for SU(2)l or 

2n3 = 24, (3.11b) 

color-triplets for SU(3)c; cf. Table 2.6. The factor 2 takes into account the corre- 
sponding conjugate multiplet as the initial state. The sum over the isospin/color 
degrees of freedom reads 

El^S-|' = ^(^«'-l) SU(iV")- (3.12) 

^This is correct since there is no factor included in the definition for the hard production rate (3.1), 
which counts the internal degrees of freedom. 



>ab 



(3.9) 
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For U(1)y, we have T^^^jj = 5ij5aiY^ / A with Yi given in Table 2.6. All matter 
fields of the MSSM carry non-vanishing hypercharge. Thus, the correct multipli- 
city factor for the U(1)y contributions is obtained by summing over all squared 
hypercharges of either scalars or fermions x\ listed in Table 2.6. We find 

2ni = 2 -^ = 22 , (3.13) 

where the factor of 2 accounts for the corresponding antiparticles. 

Processes I, J Since I and J describe the same physical processes after replacing the 
incoming fields by their conjugates, the multiplicities (3.11a), (3.11b), and (3.13) 
apply without the factor of 2. 

The diagrams of Fig. 3.2 fall into three separate classes depending on the number of 
bosons and fermions involved in the inital and final state. The BBF processes A, C, 
and J have two bosons in the initial state and one fermion in addition to the gravitino 
in the final state. Correspondingly, B, D, E, and H are BFB processes and F, G, and I 
are FFF processes. The hard part of the thermal gravitino production rate for the full 
standard model gauge group SU(3)c x SU(2)l x U(1)y can then be written as 



1 f di}p 



hard 2(27r)3£;7 4n 

3 

X 

a=l 



{2t:)H\Pi + P2-P-P:,) (3.14) 



1^1 {2iTy^2Ei 

(/bfbIM^^p^bP + /bbfIM^'^b^pP + /pFF|M(^)p|2) G(|pi - psi - r 

a=l 

with the shorthand notation 



/bfb = fB{Ei)ME2) [1 + /B(i^3)] , (3.15a) 
/bBF = fB{Ei)fB{E2) [1 - fFiEs)] , (3.15b) 
/fff = h{Ei)h{E2) [1 - /F(i^3)] • (3.15c) 
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Contributions from SU(2)l («= 2) and SU(3)c («= 3) 

For the weighted sums, we find 



(") 1 2 



BFBl 



IM, 



(a) 1 2 

bbfI 



(a) 1 2 



FFFI 




(3.16c) 



Equation (3.16c) can be rewritten as 

Ml\Agl{Nl-l) 



{") 1 2 



FFFI 



1 + 



3m?- 
G / 



Ml 



2s \ 

t s + t 



t ~^ s + t 



Ha 

2 



(3.17) 



Since s + t + u = ^rnf, where rrii are the masses of the external particles, and 
u = {Pi — -P)^, we can write s + 2t = t — u in the high-energy limit, T ^ rrii. Therefore, 

^2 



t s + t 



t u 



(3.18) 



Since the difference t — u and 1/t — 1/ii is odd under exchange of Pi and P2 , the integral 
over such terms will be zero, as long as the remaining integrand and the measure is 
even under this transformation. Therefore, there is no contribution to the integral from 
s + 2t in (3.16b) and we can further substitute s by —2t in the last term of (3.16b). In 
Eq. (3.17) we use the following replacements 



S2 




S^ 




+ 




t 


s + t 


S^ 




s^ 




+ 




t 




s + t 



0, 



2£2 
t 



(3.19a) 
(3.19b) 



For the gauge groups SU(2)l {a = 2) and SU(3)c (a = 3), we can therefore write the 
weighted sums of the squared matrix elements of Table 3.2 in terms of three distinct 
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matrix elements, 



(") 1 2 



BFBl 



(a) 1 2 



BBFI 



(") 1 2 



FFFI 



1 + 



1 + 



1 + 



Ml \ 


J / n -rO 


1) 


3mgy 






Ml\ 


'^aliNl- 


1) 


3m|y 








^gliNl- 


1) 


3m|; 





I Ml 



" 2 



(3.20a) 
(3.20b) 
(3.20c) 



with defined in (3.11) and 



iMil 



I Mo 



|M,M = - . 



2s 



2s' 



(3.21a) 
(3.21b) 
(3.21c) 



Recall from (3.10) and (3.12) that iV^ denotes the dimension of the fundamental repre- 
sentation of the corresponding gauge group, i.e., N2 = 2 and = 3. 

Contribution from U(1)y {a= 1) 

In analogy to the considerations which lead to (3.20), we obtain for the contributions 
from the U(1)y interactions 



IM 



(1) |2 



BFBl 



IM 



(1) |2 



BBFl 



IM 



(1) |2 



FFFI 



1 + 



1 + 



1 + 



M 



,t2 



3m~ 

G / 



M|1(|M,P-2|M,P) 



3m~ 

G / 

G / 



Ml 

4gini 
Ml 



IM, 



|Mi| 



2IM2P) 
IM3P) , 



(3.22a) 
(3.22b) 
(3.22c) 



with ni as defined in (3.13) and gi = g' . Note that Mi in the prefactor of (3.22) denotes 
the Bino mass while the squared matrix elements |Mjp are the ones of (3.21). 

For the calculation of the gravitino production rate (3.14), various integrations are nec- 
essary. Since the structure of the integrands is identical for all factors Qa of the standard 
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Table 3.2: Multiplicity coefficients for the squared matrix elements in (3.23) 



'^BFB,i a — i 

i = 1 4ni 

i = 2 -8ni 
i = 3 



a = 2 
4(iV| - 1) {N2 + f ) 
-4(7V| - l)n2 




4(iV|-l) (iV3 + ^) 
-4(iV| - l)n3 



model gauge group, we introduce the following notation 



d^G 



(BFB) 

hard 
(FFF) 

hard 
(BBF) 

hard 



a=l 
3 

3 



•^FFF.i -'fFF 



G/ 1=1 
3 



3m?- , 
G/ 1=1 



/ ^ I "^BBF.i -'bBF 



(3.23a) 
(3.23b) 
(3.23c) 



The coefficients c 



,(°) 

BFB.i 



are given in Table 3.2. The other factors are given by 



J") _ J") _ 

•-BBF, 1 ~ ''FFF, 2 ~ 

J") _ J") 
'-BBF, 2 ~ ''BFB, 2 ' 

_ _ 
''BBF, 3 ~ ''FFF, 1 ~ 



-BFB, 3 ' 



J") _ J") 
"''FFF, 3 ~ ''BFB,1 



(3.24a) 
(3.24b) 
(3.24c) 



The calculation of the integrals , I^^p , and /g^gp is presented in detail in Ap- 

pendix C. Here are the results: 



r|A/iP 
' r BFB \ 
I FFF J 



+ 



1 ] T^ME) 
1/2 J 1927r4 

1 



k* 6 



7 + 



C'(2) 



In 



C(2) 
El — E3 
E3 



In 2 





E 



rBFBT. + £^1) 

I FFF J ' 



+ e(^ - El) 



fs BFB 1 
I FFF J 



MM 



Ei 



(3.25) 
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xleiE- i?3)|| - E^) + e{E, -E)(^- E2 



E"^ \3 J ^ ^ V 3 

+ 9(^3 - E2)e{E - E3) ^\~f'' [{E2 - Es){E2 + 2^3) - 3(^2 + ^3)^] 

- G{E2 - E)e{E - Es) ^^'~^P' 

{E2 - Ef 



+ Q{E - E2)e{E3 - E) 



3^2 



Q{E2 - Es)e{E3 - E)^^^^ [{E2 - Es){E2 + 2^3) - 3(^2 + ^3)^] 



(3.26) 



^ noo rE+Ea 



x{Q{E- Es) 



1 EfE? 



3 



E"^ E + E3 



+ e{Es-E)- ^2 



E + £^3 

- QiE - E3)Q{Es - E2)^^^^ [£2(^3 -E)- E^iE^ + E)] 

+ Q{E3 - E)Q{E2 - £3) "^'^2^' [^2(^3 -E)- Es{E^s + £)] | , (3.27) 

where 7 = 0.57722 ... is the Euler-Mascheroni constant, C{x) is Riemann's Zeta function 
with C(2) = 7r2/6 and C'(2)/C(2) = -0.56996 ... . In the above expressions, Ei + E2 = 
E3 + E is understood. 

Only the part oc 1/t in the squared matrix element |Mi|2 given in (3.21a) exhibits a 
singular behavior for t ^ 0. Therefore, the integrals contributing to the logarithmic 
dependence on k* are II \ • The logarithm is extracted analytically using integration 

I FPF J 

by parts. Details are given in Appendix C; see Eq. (C.31) and subsequent steps. In the 
calculation of /I^bfb \ ^iid /l^Fp , the cutoff k* is set to zero from the very beginning. 

I BBF J I BBF J 

The hard part of the gravitino production rate (3.14) is then given by the sum 

(BBF) 

hard '^^P hard 



hard 



+ ^ 



(^^^) dV^ 



d^p 



(3.28) 

hard 
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We correct an error in the SU(3)c result of Ref. [2] 



We do not find the term 



r3(iV + nj)[Li2(-e-^/^) - vrVe] 



(3.29) 



given as part of Ibfb in (C.14) of Ref. [2]. Although such a contribution appears as 
a surface term in our calculation, it is canceled by another surface term. The crucial 
spots to look at are Eqs. (C.42) and (C.49) in our Appendix C. The authors of [2] agree 
with our finding and will publish an erratum. The phenomenological implications will 
be discussed in chapters 4 and 5. 



As mentioned in Sec. 3.1, hard thermal loop (HTL) resummed propagators have to be 
used for soft gauge boson momentum transfers. Since diagrams of nominally higher 
order in the loop expansions can contribute to same order in the coupling constant g at 
high temperatures, this corresponds to an improved perturbation theory where a certain 
subset of diagrams, the HTL self-energies, are resummed [19]. 

The HTL self-energies which contribute to the gauge boson polarization tensor n|j"^ are 
given by one-loop diagrams such as the ones shown in Fig. 3.1. They carry soft external 
momenta K = (A:o,k) with ko and A; = |k| of order gT and hard internal loop momenta 
of order T. In the HTL approximation of hard loop momenta, only a small part of the 
integration region in the one-loop diagrams contributes. Hard thermal loops are gauge 
invariant and the corresponding self-energy satisfies the Ward identity K^H^f^ (K) = 0. 
They are exclusively due to thermal fluctuations and are therefore ultraviolet finite [16]. 

At finite temperatures, the thermal bath constitutes a privileged rest frame. The polar- 
ization tensor H^j"^ can be decomposed into two independent propagating modes which 
are chosen to be the longitudinal and the transverse parts of the self-energy of the gauge 
bosons Both are physical and read respectively 



In the high temperature limit, kQ,k <^ T, to leading order in the gauge couplings ga, 

^For a direct comparison with Ref. [2], it should be stressed that our definition for the production 
rate differs by a factor of /p(ii')/(27r)^. 



3.3 Soft Contribution 



(3.30a) 




(3.30b) 
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these components of the self-energy are given by [20] 

lt\k,,k) = -?,ml (l-|ln^l , (3.31a) 



1 _ ( 1 _ _ ] ^ In ^0 + ^ 
kl) 2k ^ ko- k 



(3.31b) 



where denotes the thermal masses of the corresponding gauge boson A^"^ . Note that 
the self-energy components (3.31) depend on the gauge boson momenta in a non-trivial 
way and that they have an imaginary part for k^ < k"^. 

In the static limit, /cq 0, the longitudinal part of the self-energy reduces to n[^°^^(0, k) = 
— 3m^ which is analogous to the Debye screening of static electric fields in a QED plasma 
with inverse screening length A^^ = VSma- There is no static magnetic screening since 
Ili^\o,k) = 0. Nevertheless, the transverse self-energy IIt approaches zero sufficiently 
slow so that quantities in which the magnetic divergence is only logarithmic are regu- 
larized [14]. 

In SUSY extensions of the Standard Model, the hard thermal loops shown in Fig. 3.1 
contribute to the self-energy tensor and thereby to the thermal gauge boson mass nia, 

Tfl^ gauge bosons ~^ gauginos ~^ '^a, fcrmions ~^ '^cijScalars • (3.32) 

For the non-abelian gauge group SU(iVQ), the contributions are given by [16, 21] 

"^a, gauge bosons — g 1/^ ^ I = ^ , (3.33a) 

b,c 



2 ya 
m — 
""q, gaugmos 



2rp2 2rp2 

b,c 



18 ' 18 



"^a, fermions ~ ^ 



2 



EKl\' = na'-^, (3.33c) 



18 ''•^J' 36 



_ 9a^ \ ^ irrja) |2 _ 9a^ 

^a, scalars 



2 



For a = 2 and a = 3 with defined in (3.11), this determines the thermal masses of 
the wino and the gluino, respectively. They read 

< = '-^(No. + ^). (3.34) 



6 V " ■ 2 

Analogously, one finds for U(1)y: 



2 afniT^ 

m| = — !— (3.35) 
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(a) Gauge boson contributions for SU(3)c and SU(2)l 




(b) Gaugino contributions 




(c) Fermion contributions 



K 



K 



(d) Scalar contributions 



(a) 



Figure 3.2: HTL self-energy contributions to the gauge boson polarization tensor Yv^J 
The internal loop momenta are hard, i.e., of order T. In the resummation these contri- 
butions become relevant for soft gauge boson momenta K ~ gT. 



which determines the thermal bino mass. In the U(1)y case, m\ gauge bosons ~ ^ ^''■^^ 



m 



1 , gauginos 



since the diagrams of Fig. 3.1a are absent. 



In covariant gauge, the HTL-resummed gauge boson propagator has the form [22, 23] 



with the tensorial quantities 
A 



2 



K * V ( K • V 



2 • 



(3.36) 

(3.37a) 
(3.37b) 
(3.37c) 
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The velocity of the thermal bath is denoted by v and the gauge fixing parameter by ^. 
The transverse and longitudinal propagators are 

A^"^(A;o, k) = '^—r-. , (3.38a) 

^ kl-k^-I['^^\k^,k) 



which have spectral representations [22] 

du:^-—p^^]^{uj,k) . (3.39) 
For \uj\ < k, the spectral densities p^^^rp are given by 

Pt^ i^, k) = (1 _ ^2) [^^(^)2 + + Bt{xW] ' ^^■^^''^ 
( ) 3 2x 

with X = uj/k and z = k"^ /m^ and 



^t(x) = ^TTX , i?T(^) = 7 f2-^ +xlnl±^ ) , (3.41a) 

4 4\1 — 1 — X 



Ai^ix) = ^TTX , Bi^ix) = ^(^2 - xlni±|) . (3.41b) 

Let us now turn to the calculation of the soft part of the gravitino production rate. 
In the previous section we have seen that the squared matrix elements for the 2 — > 2 
scatterings given in Table 3.2 contain the factor 

|MPoc(l + i|). (3.42) 

The first term results from the helicity ±3/2 states of the gravitino, the second term 
from the helicity ±1/2 states of the gravitino which represent the goldstino components. 
For the gravitino self-energy Eg, it has been shown up to two loop order in the gauge 
couplings that one obtains the same factor (3.42); cf. [2]. Indeed, we employ the effective 
theory for light gravitinos given in (2.83) to calculate the production rate for the helicity 
±1/2 components of the gravitino in terms of the imaginary part of the goldstino self- 
energy s|f*^(P), namely. 



dr 



(a) 



soft 



—fp{E)lm^"'{E + ie,p)\k<k' . (3.43) 
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A(") 




Figure 3.3: The leading contribution to the imaginary part of the goldstino self-energy 
in the effective theory of light gravitinos. The blob indicates the HTL-resummed gauge 
boson propagator (3.36). 

This leads to the full rate by replacing M^/3m~ with the prefactor (3.42). Note that we 
have introduced the cutoff k* since we consider here only soft three- momentum transfers. 

The leading contribution to the self-energy is given by the gauge boson-gaugino loop 
shown in Fig. 3.2; see the discussion below (2.83). Here the HTL-resummed gauge 
boson propagator (3.36) is indicated by the blob. Accordingly, one finds for the leading 
contribution to the self-energy^ 

-'^f'(^) = {^'c-^ V^-n I'^fc' ir-f] . 

(3.44) 

Working in the high-energy limit, we have neglected the masses in the gaugino propaga- 
tors. The polarization sum for the goldstino spinors u^{P) and u^{P) is included. The 
gauge bosons carry soft momenta K and the gaugino momenta are given hy Q = P — K. 
The multiplicity factors X^"^ count the number of gauge bosons/gauginos in the loop 
for the corresponding gauge group Qa, namely, 

= iV^ - 1 for SU(iV„) , (3.45a) 
X(^) = 1 for U(1)y . (3.45b) 

One computes (3.44) in the imaginary time formalism where the energies of the particles 
are given by their Matsubara frequencies. For the gauge bosons, the frequencies are 

^We commit a certain abuse of terminology since a self-energy contribution is usually referred to as 
the expression which one obtains from the amputated diagram of Fig. 3.2. The inclusion of the spinors 
u{P) and u{P) is necessary to interpret (3.43) in terms of a probability (see [15]). 
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even, ko = 2mnT, and the integral over ko turns into a sum over discrete energies, i.e., 

^-■^ E ■ (3-40) 

n=— OO 

After performing the polarization sum, the self-energy reads 

d^k 1 



(,) 4x(^)<r^ 



(27r)3 g2 



with the Dirac traces: 



1 



tr{i^[^^7l ^[^^7^]V} > 



(3.47) 

(3.48a) 
(3.48b) 



Note that the dependence on the gauge fixing parameter drops out since (3.37c) 
contracted with the gauge boson momentum K vanishes in the Dirac trace. 

Using the spectral representations of the propagators (3.39), the summation over the 
Matsubara frequencies can be performed conveniently with the Saclay method [24] . The 
result reads after analytic continuation from discrete energy values po to continuous real 
goldstino energies E [2] 



(a) 



d^p 



soft 



G 



h{E) 



k* 



dkk-n - 

-k ^ 



Pl(w, k)[l--g^]+ pt(^, ^0(1- — 



k^ 



(3.49) 



The structure of the integrand is identical to the one obtained for the axion produc- 
tion rate [14] where the analytic dependence on the cutoff k* is extracted. Thus, after 
performing the integrations, one finds 



dr 



(a) 



d'^p 



soft 



^^^^^ WM>r 

G 



In 



1.379 



(3.50) 



result for the soft part of the gravitino production rate (3.6) 



The replacement of M^/3m~ with the factor (3.42) for the full theory yields the final 

G 



d^p 



soft 



a=l 



3m1 

G / 



327r4M2 



In 



*2 



rat. 



1.379 



(3.51) 



with the multiplicites X^'^'^ given in (3.45) and the thermal masses given in (3.34) 
and (3.35). 

Adding the results for the hard and the soft part of the gravitino production rate, one 
finds that the logarithmic dependence on k* cancels out. We will do so by calculating 
the Boltzmann collision term which is the crucial quantity for all further calculations. 
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3.4 The Boltzmann Collision Term 



The quantity we are interested in is the gravitino number density'^ 



ng = 4 



(2^)3 

Its evolution with cosmic time t is governed by the Boltzmann equation 

dnp 



+ 3Hn^ = . 
dt G G 



(3.52) 



(3.53) 



The second term on the left-hand side accounts for the dilution of gravitinos due to 
the expansion of the Universe, which is described by the Hubble parameter H. For 
negligible gravitino disappearance processes, the collision term Cg on the right-hand 
side of the Boltzmann equation describes gravitino production processes. It is obtained 
by integrating the thermal gravitino production rate 



d^p 



d^p 



(fp 



d^p 



hard ^^P 



soft. 



(3.54) 



Consider first the soft part (3.51). The corresponding collision term reads 



G, soft 



1 + 



M2 \ 3X(°)m2 r 



3m?- 
G / 



In 



„*2 



mi: 



1.379 



dE , — - 

e^/T + 1 



(3.55) 



The final integration over the Fermi-Dirac distribution function is easily performed. The 
result reads 



^G,soft 



G/ 



167r3M2 



In 



1.379 



(3.56) 



Now we turn to the hard part (3.28). The production rates for the BFB, FFF, and BBF 
processes are given in (3.23). The numerical integrations are performed using VEGAS [25]. 
We get: 



^BFB 

G.hard 



E 



9lT' L , Ml 



3?Tl~ 

G , 



» 

-BFB, 1 



|^(,n(?^ I +0.9930 



11.1362 X 10"^ 



^ '^BFB,2 



1.3284 X 10" 



(3.57) 



■^The factor of four accounts for the internal degrees of freedom of the massive gravitino. 
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"-"Chard Mi 



1 + 



Ml 

3m?- 
G / 



X < c 



(a) 

FFF,1 



6.9992 X 10" 



+ c- 



(a) 

EPF, 3 



0.5039 X 10-^ 



(3.58) 



^BBF 

G, hard 



2^6 



Ml 



(a) 
^BBF, 2 



1 + 



3m?- 
G / 



1.2975 X 10" 



+ c- 



(a) 

BBF,3 



0.8647 X 10" 



(3.59) 



In the sum of the soft and hard parts, the logarithmic dependence on k* cancels out as 
anticipated. Plugging in the squared thermal masses (3.34) and (3.35), the final result 
can be brought into the following form 



C. 



G 



+3m| j levrSMl"^"^" U« 



Q = l 



G, 



(3.60) 



This is one of the main results of this thesis. The coefficients are given by Cq, = 
(11, 27,72) and the scales in the logarithms by ka = (1.266, 1.312, 1.271). These 
numbers are associated with the gauge groups U(1)y, SU(2)l, and SU(3)c, respectively. 
The temperature T provides the scale for the evaluation of the gaugino mass parameters 
Ma = {Ml, M2, M3 ) and the gauge couphngs ga = (g' , 9 , Os)- 

The error in the hard production rate of Ref. [2] manifests itself in a larger coefficient 
k^ = 1.271 for the SU(3)c contribution than 1.163 obtained from [2]. 

Recall that the Braaten-Yuan prescription [14] relies on the weak coupling limit g <^ 1. 
Because of the large value of the strong coupling constant, e.g., (7s(100GeV) ~ 2.5, the 
results for the production rates and the corresponding collision term require high temper- 
atures T » 10^ GeV, where, for example, gs{10^ GeV) ~ 0.99 and g^iW^ GeV) ~ 0.88. 

Since we neglect gravitino disappearance processes, the collision term (3.60) acts as pure 
source term in the Boltzmann equation (3.53). For small temperatures, the logarithm 
turns negative for the SU(3)c part. Thus, the result becomes unphysical and shall not 
be trusted when extrapolating to small temperatures. 
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Chapter 4 



Gravitino Cosmology 



The new result for the Boltzmann colhsion term (3.60) has important imphcations for 
gravitino cosmology. In this chapter, we compute the gravitino yield which describes 
the primordial gravitino abundance. This quantity is crucial for phenomenological con- 
siderations of stable and unstable gravitinos. For gravitino dark matter scenarios, we 
calculate the relic gravitino density from thermal production. This consequently yields 
an upper bound on the reheating temperature of the Universe. 

4.1 Gravitino Yield from Thermal Production 

The observed isotropy and homogeneity of the Universe on large scales allows us to 
express the overall geometry of the Universe in terms of the Robertson- Walker metric 
with the line element 



where {t,r,6,(j)) are comoving coordinates. By a rescaling of r, the curvature param- 
eter k can be assigned the discrete values k = 1, —1, or 0, corresponding to spatially 
closed, open, or flat geometries. The evolution of the scale factor R is described by the 
Priedmann equation 



which defines the Hubble parameter H. The total energy density of the Universe is 
denoted by p. The derivative of the scale factor R with respect to cosmic time t is 




(4.1) 




(4.2) 



written as R. 
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4.1. Gravitino Yield from Thermal Production 



In the radiation dominated epoch of the Universe, p is given in good approximation by 

P = 9*%T\ (4.3) 

where T is the photon temperature and denotes the effectively massless degrees of 
freedom, i.e., those species with mass rrii <^Ti, 

i=bosons i=fcrmions 

The Hubble rate in the radiation dominated epoch is given by {k = 0) 



Hm = J^— , (4.5) 

where Mp is the reduced Planck mass (2.4). During this epoch, time and temperature 
are related via H{T) = l/{2t). 

The entropy density of the Universe, defined as s = {p + p)/T, is dominated by the 
contribution of relativistic particles for which p = p/3 holds. Hence, one finds 

s = g.s^T\ (4.6) 

where 

j=bosons i=fermions 

In the previous section, the Boltzmann equation (3.53) has been written in terms of the 
gravitino number density ng. It is useful to scale out the expansion by dividing the 
gravitino number density by the entropy density s. This defines the yield variable: 

^ ^ . (4.8) 

^ s 

With the conservation of entropy per comoving volume, sR^ = const., the Boltzmann 
equation (3.53) can be rewritten as 

— ^ = ^ . 4.9 

dt s ^ ' 

Using dt = —dTf[H{T)T], the gravitino yield is obtained by integrating 

Cr{T)dT 

= -li^ ■ ("») 
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4.2. Gravitino Dark Matter 



We consider thermal gravitino production beginning after completion of the reheating 
phase where the temperature of the primordial plasma is the reheating temperature Tr. 
We assume that any initial gravitino population has been diluted away by inflation, i.e., 
Yq(T^) = 0. Hence, the gravitino yield at the temperature T' is 



Unstable gravitinos have typically long lifetimes because their interactions are suppressed 
by Mp. In particular, when the gravitino is lighter than < 20TeV, its lifetime becomes 
longer than > Is [26]. Hence, unstable gravitinos may decay during and/or after big- 
bang nucleosynthesis where tasN —Is and Tbbn — 1 MeV. Thus, the yield of gravitinos 
from thermal production prior to their decay is obtained for T' = T^bn (it^q ~ 20TeV). 
Note that the dependence of the collision term (3.60) cancels out in the integrand 
of (4.11). Furthermore, recall that we consider scenarios in which Tr > 10^ GeV ^ Tbbn 
(see Section 3.4). Thus, with the collision term (3.60), we can solve the Boltzmann 
equation to good approximation analytically. We find 



with ya = (0.653, 1.604, 4.276) x lO'^^ foj. u(1)y, SU(2)l, and SU(3)c, respectively. 
The scales ka in the logarithms are given at the end of Sec. 3.4. Here, we have used that 
after reheating, at temperature Tr, all particles of the MSSM are in thermal equilibrium 
and relativistic, for which g^iT^i) = g*s{TR) = 915/4. 

The yield (4.12) is the starting point for studies of cosmological constraints in scenarios 
with unstable gravitinos. In the remainder of this chapter we will consider the case of a 
stable gravitino. 

4.2 Gravitino Dark Matter 

In the following, we focus on scenarios in which the gravitino is the LSP and stable due 
to R-parity conservation. Gravitinos, once produced, will thus contribute to the present 
value of the energy density p since they do not decay. Because of the large redshift, the 
gravitino energy density at the present time to is Pg(io) = i^Qi'^o)- We have seen that 
the thermal production of gravitinos is efficient only during the very early hot radiation 
dominated epoch so that 




(4.11) 



a=l 



3 




Ygin) yg(TBBN) 



(4.13) 
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4.2. Gravitino Dark Matter 



for stable gravitinos. Here, 

To = 2.725 K = 2.348 x 10"^^ GeV (4.14) 
is the present temperature of the Universe [5]. 

Thus, the present day density parameter of thermally produced gravitinos is given by: 

Pc{to) Pc{To) 

Here, the dimensionless quantity h is used to parameterize the Hubble constant Hq = 
H{Tq) = lOO/ikms""*^ Mpc^^ and the present value of the critical density reads"*^ [5] 

p^{tQ)/h^ = 8.096 X 10"^^ GeV^ . (4.16a) 

The entropy density s(To) is obtained from (4.6) with g*s{To) = 43/11. 

Thus, from (4.12), (4.13), and (4.15) we find the result for the relic density of thermally 
produced gravitinos to leading order in the Standard Model gauge couplings: 

(4.17) 

with = (0.018, 0.044, 0.117) for U(1)y, SU(2)l, and SU(3)c, respectively. 

The relic gravitino density ^l^h? is essentially linear in the reheating temperature Tr. 
Recall that the production of the helicity ±3/2 states of gravitinos is described by the 
first term in the generic factor (l + M^/3m~). Thus, the relative weights for the helicity 
±3/2 production of gravitinos for the different factors of the Standard Model gauge group 
are basically given by Ua gaiTn)'^ — which is model- independent. The production of the 
helicity ±1/2 states depends on the ratio of squared gaugino masses to gravitino 
mass niQ. Hence, the thermal production becomes more efficient for light gravitinos. 

In order to calculate numbers from (4.17), we need to evaluate the gauge couplings Qa and 
the gaugino mass parameters at the scale provided by the reheating temperature Tr. 
At the one-loop level, the renormalization group (RG) equations for the gauge couplings 
read 

'^"^"^^ - ^"%.(Q) (4.18) 



dlnQ/Qo 167r2 



^ft = 0.73tro^ [5] 
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10^ 10'' 10** 10" lO^-* Mgut 10^ lO'' 10** 10" 10" Mqut 

Q [GeV] Q [GeV] 

(a) (b) 



Figure 4.1: One-loop renormalization group running of the gauge couplings (a) and 
the gaugino mass parameters (b) in the MSSM. The point of gauge coupling unification 
is denoted as A'/qut- The runnings of the gaugino masses are inferred from (4.21) 
for universal A/1,2,3 = ™i/2 A/qut with a value of = 400 GeV. They are 

given in (b) as solid lines. Dashed lines in (b) show a non-universal scenario with 
0.5 Ml, 2 = M3 = TOi/2 at A/gut, i-c, a;i,2 = 2. 



where Q is the scale of evaluation and Qq is some input scale. One can solve (4.18) 
analytically: 



gaimzY 



In 



Ik 

mz 



-1/2 



(4.19) 



,(1) 



with Qq = mz ~ 91.19 GeV [5] and Q = Tr. The beta-function coefficients (3a 
( 11 , 1 , — 3 ) correspond to Qa = {g' , g , gs) the MSSM, respectively. In Fig. 4.0a the 
point of gauge coupling unification is referred to as the grand unification (GUT) scale 
Mgut — 2 x 10^^ GeV. It is defined as the point where the GUT-normalized hypercharge 
coupfing, 



gi 



^g' 



(4.20) 



the weak coupling, g, and the strong coupling, g^, meet. 

The one-loop RG equations for the gaugino masses Mq, are given by the expression 
analogous to (4.18) and the same coefficients (3a^ as above. Often, universal boundary 
conditions are considered for the RG equations in which the gaugino masses unify at 
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4.3. Upper Bounds on the Reheating Temperature 



Mgut- Since our result for the rehc gravitino abundance depends on Mq, we will 
also consider scenarios in which the gaugino masses do not unify at Mgut- We can 
parameterize this by writing 

1 Mi(Q) ^ 1 M2(Q) ^ MsiQ) 

XI gi{QY X2 g2{Q? 93(Q)2 ^ ■ ^ 

which holds at any scale Q.^ For the gaugino masses, we choose the input scale Qo to 
be Mgut- This defines the gaugino mass parameter mi/2) namely, 

mi/2 = M3(Mgut) = M2(Mgut)/x2 = Mi(Mgut)/xi , (4.22) 

so that for a unifying scenario, i.e., for xi^2 = 1) we have mi/2 = Mi(Mgut) = 
M2(Mgut) = M3(Mgut)- As illustrated for mx/2 = 400 GeV, fixing mii2 determines 
all gaugino mass parameters M^ . 

Taking the RG evolution into account, we compute the relic gravitino density. Fig- 
ure 4.1 shows Vi^h? as a function of the reheating temperature Tr for gravitino 
masses mg = 1, 10, 50, and 300 GeV. We consider two representative values of the 
gaugino mass parameter mi/2 = 400 GeV and 1500 GeV in Figs. 4.1a and 4.1b, respec- 
tively. The solid curves represent the gravitino density for universal gaugino masses 
-^1,2,3 = &t Mgut while the dashed curves show a non- universal scenario where 

0.5Mi^2 = M3 = mii2 at Mgut, i-e-, xi^2 = 2. The corrected result for the SU(3)c 
contribution is given by the dotted lines. The grey band indicates the dark matter 
density [5] 

^Arah^ = 0.105t[!;0?^ . (4.23) 

and thus shows the parameter region in which the thermally produced gravitinos provide 
the observed dark matter density. 

We find that electroweak processes enhance by about 20% for universal gaugino 
masses at Mgut- In non-universal cases, Mi^2 > M3 at Mgut, the electroweak con- 
tributions are more important. For 0.5Mi^2 = M3 at Mgut, they provide about 40% 
of ^ij^h?. Moreover, with our new value — see Eqs. (3.60) and (4.17) — we find an 
enhancement of about 30% of the SU(3)c contribution to relic density in comparison to 
the result given in [2]. 

4.3 Upper Bounds on the Reheating Temperature 

In the gravitino dark matter scenario we can derive an upper bound on the reheating 
temperature Tr from < fidm, once mi/2 is specified. Figure 4.2 shows the upper 

^Up to small two loop effects and possible (unknown) threshold effects close to Mgut [9]. 
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4.3. Upper Bounds on the Reheating Temperature 



limits on Tr for mx/2 = 400 GeV and mi/2 = 1500 GeV, respectively. The solid and 
dashed curves give the upper bounds on Tr inferred from our SU(3)c x SU(2)l x U(1)y 
result of Qj^h"^ for universal {Mi 2.3 = 1^^1/2) ^-i^d non-universal (0.5 Mi 2 = M3 = 
gaugino masses at Mgut? respectively. The dotted curves show the SU(3)c limits for 
M3 = mi/2 at Mgut- We have adopted 

^diT^^ = 0.126 (4.24) 

as a nominal 3cj upper limit on fidm^^- 

Note that for higher values of the bounds on Tr are more stringent; cf. (4.17). For 
small mg, the thermal gravitino production is very efficient. Then J^^^/i^ ~ M^/3mQ 
since the production of helicity ±1/2 states dominates. For large values of mg, the 
upper limit on Tr becomes more severe again. Then Q.J^h^ ~ mp; since the production 
of the helicity ±3/2 states constitutes the dominant part. The bounds on Tr will become 
more severe when one includes non-thermal production processes such as late-decays of 
the next-to-lightest supersymmetric particle (NLSP) into the gravitino. This, however, 
depends on the details of the realized SUSY model while the bounds in Fig. 4.2 are 
rather model independent. 
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4.3. Upper Bounds on the Reheating Temperature 




Figure 4.2: The relic gravitino density from tliermal production, flj^h^, as a function 
of Tr. The solid and dashed curves show the SU(3)c x SU(2)l x U(1)y results for 
universal (Mi_2,3 = fni/2) and non-universal (0.5 Mi_2 — -^^3 = ^1/2) gaugino masses at 
AfcuT, respectively. The dotted curves show the new result of the SU(3)c contribution 
for A/3 = 7711/2 at A/gut- The grey band indicates the dark matter density fJdm^^- 
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Figure 4.3: Upper bounds on Tr from Q~ < The solid and dashed curves show 

the Umits inferred from the fuU SU(3)c x SU(2)l x U(1)y result (4.17) for Mi,2.3 = mi/2 
and 0.5A/i.2 = — mi/2 at A/qut, respectively. The dotted curves show the limit 
on Tji for our new result of the SU(3)c contribution. 
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Chapter 5 

Testing Leptogenesis at Colliders 



The smallness of the neutrino masses can be understood naturally in terms of the see-saw 
mechanism [27, 28] once the Standard Model is extended with right-handed neutrinos 
which have heavy Major ana masses and only Yukawa couplings. For a reheating tem- 
perature after inflation, Tr, which is larger or not much smaller than the masses of the 
heavy neutrinos, these particles are produced in thermal reactions in the early Universe. 
The CP-violating out-of-equilibrium decays of the heavy neutrinos generate a lepton 
asymmetry that is converted into a baryon asymmetry by sphaleron processes [29] . This 
mechanism, known as thermal leptogenesis, can explain the cosmic baryon asymmetry 
for Tr > 3 X lO'^GeV [30]. 

One will face severe cosmological constraints on Tr if supersymmetry is discovered. 
We have seen that gravitinos are produced efficiently in the hot primordial plasma. 
Because of their extremely weak interactions, unstable gravitinos with mg < 5 TeV 
have long lifetimes, rg > 100 s, and decay after BBN. The associated decay products 
affect the abundances of the primordial light elements. Demanding that the successful 
BBN predictions are preserved, bounds on the abundance of gravitinos before their decay 
can be derived which imply Tr < 10^ GeV for rriQ < 5 TeV [26]. Thus, the temperatures 
needed for thermal leptogenesis are excluded. 

Let us therefore consider SUSY scenarios in which a gravitino with mg > 10 GeV is the 
LSP and stable due to R-parity conservation. These scenarios are particularly attractive 
for two reasons: (i) the gravitino LSP can be dark matter and (ii) thermal leptogenesis 
can still be a viable explanation of the baryon asymmetry [31]. 
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5.1. Collider Predictions of Leptogenesis 




mg [GeV] 

Figure 5.1: Upper limits on the gaugino mass parameter TOi/2 from < i^™,^^ for 
the indicated values of Tr. The solid and dashed curves show our SU(3)c x SU(2)l x 
U(1)y results for universal (A/1.2,3 = ^1/2) non-universal (0.5Afi,2 = M3 = m.1/2) 
gaugino masses at il/cuT, respectively. The dotted curves show the SU(3)c limits for 
M3 = mi/2 at Mgut- 



5.1 Collider Predictions of Leptogenesis 

Thermal leptogenesis requires Tr > 3 x 10^ GeV [30]. This condition together with the 
constraint fi-^ < ^^^^ [see Eqn. (4.24)] leads to upper limits on the gaugino masses. 
The SU(3)c result for fl^ implies limits on the gluino mass [2, 32]. With our SU(3)c x 
SU(2)l X U(1)y result, the limits on the gluino mass M3 become more stringent because 
of the new value (3.60) and the additional electroweak contributions. Moreover, as 
a prediction of thermal leptogenesis, we obtain upper limits on the electroweak gaugino 
mass parameters Mi^2- At the Large Hadron Collider (LHC) and the International 
Linear Collider (ILC), these limits will be probed in measurements of the masses of the 
neutralinos and charginos, which are typically lighter than the gluino. If the superparticle 
spectrum does not respect these bounds, one will be able to exclude standard thermal 
leptogenesis. 

Figure 5.0 shows the gaugino mass bounds for Tr = 10^, 3 x 10^, and 10^^ GeV evolved 
to Mgut 5 i-e., in terms of limits on the gaugino mass parameter mi/2- With the observed 
superparticle spectrum, one will be able to evaluate the gaugino mass parameters Mi, 2,3 
at Mgut using the SUSY renormalization group equations [33, 34, 35, 36]. While the 
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5.2. Decays of the Next-to-Lightest Supersymmetric Particle 



determination of Mi^2 at low energies depends on details of the SUSY model that will be 
probed at colliders [9], the bounds shown in Fig. 5.0 depend mainly on the Mj relation 
at Mgut- This is illustrated by the solid and dashed curves obtained with Mi^2,3 = "^112 
and 0.5Mi^2 = Afs = "ii/25 respectively. The dotted curves represent the SU(3)c limits 
for M3 = at Mgut and emphasize the importance of the electroweak contributions. 

5.2 Decays of the Next-to-Lightest Supersymmetric Par- 
ticle 

With a gravitino LSP of > 10 GeV, the next-to-lightest SUSY particle (NLSP) has 
a long lifetime of tnlsp ^ 10^ s [37, 38]. After decoupling from the primordial plasma, 
each NLSP decays into one gravitino LSP and Standard Model particles. The resulting 
relic density of these non-thermally produced gravitinos is given by 

^f'h' = n^LSPh' , (5.1) 

where wnlsp is the mass of the NLSP and f^NLSp/^^ is the relic density that the NLSP 
would have today, if it had not decayed. As shown below, more severe limits on mi/2 
are obtained with fl^^h^ taken into account. Moreover, since the NLSP decays take 
place after BBN, the emitted Standard Model particles can affect the abundance of the 
primordial light elements. Successful BBN predictions thus imply bounds on and 
"^NLSP [37, 38]. From these cosmological constraints it has been found that thermal 
leptogenesis remains viable only in the cases of a charged slepton NLSP or a sneutrino 
NLSP [32, 39]. 

5.3 CoUider Tests of Leptogenesis 

Thermal leptogenesis will predict a lower bound on the gravitino mass once the 
masses of the Standard Model superpartners are known. With a charged slepton as the 
lightest Standard Model superpartner, it could even be possible to identify the gravitino 
as the LSP and to measure its mass itIq at future colliders [40, 41, 42, 43]. Confronting 
the measured with the predicted lower bound will then allow us to decide about the 
viability of thermal leptogenesis. 

To be specific, let us assume that the analysis of the observed spectrum [35, 36] will point 
to the universality of the soft SUSY breaking parameters at Mgut and, in particular, 
to the minimal supergravity (mSUGRA) scenario with the gaugino mass parameter 
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^^1/2 = 400 GeV, the scalar mass parameter tjiq = 150 GeV, the trilinear coupling 
Aq = —150, a positive higgsino mass parameter, ^ > 0, and the mixing angle tan/3 = 30 
in the Higgs sector. A striking feature of the spectrum will then be the appearance of the 
lighter stau ti with m^^ = 143.4 GeV as the lightest Standard Model superpartner [44] . 
In the considered gravitino LSP case, 10 GeV < < m^f^ , this stau is the NLSP and 
decays with a lifetime of t:^^ > 10^ s into the gravitino. For the identified mSUGRA 
scenario and the considered reheating temperatures, the cosmological abundance of the 
Ti NLSP prior to decay can be computed from Onlsp^^ = ^r-^h? ~ 3.83 x 10"'^, which 
is provided by the computer program micrOMEGAs [45]. For given mg, this abundance 
determines Q^^h? and the release of electromagnetic (EM) and hadronic energy in ti 
NLSP decays governing the cosmological constraints [37, 38]. 

Figure 5.1 allows us to probe the viability of thermal leptogenesis in the considered 
mSUGRA scenario. 1 From the constraint + Og'^^ < a"^^, we obtain the solid 
curves which provide the upper limits on mx/2 for Tr = 10^, 3 x 10^, and 10^^ GeV. The 
dashed line indicates the mi 12 value of the considered scenario. The vertical solid line 
is given by m^f^ = 143.4 GeV which limits from above. In the considered scenario, 
the mi/2 value exceeds the mi 12 limits for Tr > 10^*^ GeV. Thus, temperatures above 
lO^'' GeV can be excluded. Temperatures above 3 x 10^ GeV and 10^ GeV remain 
allowed for mg values indicated by the dark-shaded (dark-green) and medium-shaded 
(light-green) regions, respectively. The values indicated by the light-shaded (grey) 
region are excluded by BBN constraints for late ri NLSP decays.^ 

Here thermal leptogenesis, Tr > 3 x 10^ GeV, predicts > 130 GeV and thus a ti 
lifetime of t^^ > 10^^ s [37, 38]. If decays of long-lived ri's can be analyzed at colliders 
giving evidence for the gravitino LSP [40, 41, 42, 43], there will be the possibility to 
determine in the laboratory: From a measurement of the lifetime t^-^ governed by 
the decay ri — > Gr, m^ can be extracted using the supergravity prediction for the 
associated partial width, 

r~ ^ t-\ti - Gt) = ^ 1 - ^ (5-2) 

m~ \ m~ ] 

^Thermal leptogenesis requires right-handed neutrinos and thus an extended mSUGRA scenario. This 
could manifest itself in the masses of the third generation sleptons [33]. Since the effects are typically 
small, we leave a systematic investigation of extended scenarios for future work. 

^We use the conservative BBN bounds considered in [38]. The average EM energy release in one ri 
NLSP decay is assumed to be Et/2, where Et is the energy of the tau emitted in the dominant 2-body 
decay t\ Gt (cf. Fig. 16 of Ref. [38]). With an EM energy release below Et/2, the grey band can 
become smaller. For less conservative BBN constraints and/or enhanced EM energy release, the excluded 
region becomes larger. 
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Figure 5.2: Probing the viability of thermal leptogenesis. The solid curves show the 
limits on the gaugino mass parameter mx/2 from + fl^^ < JIJ^^'^ for Tr = 10^, 
3 X 10^, and 10^" GeV. The dashed line indicates the mi/2 value of the considered 
scenario. The vertical solid line is given by the ri NLSP mass which limits the gravitino 
LSP mass from above: < m^:^ = 143.4 GeV. The mg values at which temperatures 
above 3 x 10^ GeV and 10^ GcV remain allowed are indicated by the dark-shaded (dark- 
green) and medium-shaded (light-green) regions, respectively. The values within 
the light-shaded (grey) region arc excluded by BBN constraints. 

Moreover, for > 0.1 m^^, can be infered kinematically from the energy of the 
tau, Er, emitted in the 2-body decay ti — > Gt [40, 42]: 

= m~ — — 2mf^ Er ■ (5-3) 

While within the dark-shaded (dark-green) region will favor thermal leptogene- 
sis, any outside of the medium-shaded (light-green) region will require either non- 
standard mechanisms lowering the Tr value needed for thermal leptogenesis or an alter- 
native explanation of the cosmic baryon asymmetry. 
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Chapter 6 



Conclusions 



The starting point of this thesis has been a discussion of the general supergravity La- 
grangian in four spacetime dimensions. We have shown how supersymmetry can be 
broken spontaneously. For a simple gravity-mediation scenario, we have explicitly car- 
ried out the transition to a softly-broken supersymmetric theory as the low-energy limit 
of supergravity. Relating the obtained effective theory with the MSSM in the high-energy 
limit of unbroken electroweak symmetry and identifying the relevant gravitino interac- 
tions, we have derived all Feynman rules necessary for the calculation of the regeneration 
of gravitinos after inflation. 

We have used the Braaten-Yuan [14] prescription together with the hard thermal loop 
resummation technique [16] to calculate the thermal gravitino production rate. We have 
considered the regeneration of gravitinos that starts after completion of reheating. As 
one of the main results of this thesis, we present the Boltzmann collision term (3.60) to 
leading order in the SU(3)c x SU(2)l x U(1)y couplings. It acts as the source term in the 
Boltzmann equation which governs the time evolution of the gravitino number density 
in the thermal bath. The collision term has been obtained in a gauge invariant way 
within the framework of thermal field theory and does not depend on arbitrary cutoffs. 
Our result includes for the first time the SU(2)l x U(1)y sector. Moreover, we correct 
an error in the previously known SU(3)c result which was obtained in Ref. [2]. 

With direct implications for gravitino cosmology, we obtain the gravitino yield (4.12) 
from thermal production by solving the Boltzmann equation. The yield parametrizes the 
primordial abundance of gravitinos in the early Universe. Hence, it is a crucial quantity 
for both scenarios with stable gravitinos and scenarios with unstable gravitinos. Focusing 
on gravitino dark matter scenarios, we have obtained the relic abundance of thermally 
produced gravitinos, ^2-^. It depends on the reheating temperature Tr, the gravitino 
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mass m^, and the gaugino masses Mq,. We find that electroweak processes enhance 
by about 20% for universal gaugino masses at the grand unification scale Mqut- For non- 
universal scenarios, the electroweak contributions become more important. Furthermore, 
with our new SU(3)c result, we find an enhancement of the SU(3)c contribution to O-^ 
by about 30% as compared to the relic density obtained in [2]. We also give an update 
for the upper bounds on the reheating temperature depending on the specified gaugino 
mass parameter mi/2- 

With the result for the relic density of thermally produced gravitino LSPs, new gravitino 
and gaugino mass bounds emerge as a prediction of thermal leptogenesis which requires 
^ 3 X 10^ GeV. If supersymmetry is realized in nature, these bounds will be accessible 
at the LHC and the ILC. For certain SUSY parameter regions, a charged slepton will 
be the NLSP. We have studied a scenario in which the lighter stau is the NLSP. Here we 
have also taken into account the non-thermal gravitino production from the stau decays 
into the gravitino LSP. There exists then the exciting possibility to identify the gravitino 
as the LSP and to measure its mass. Confronting the measured gravitino mass with the 
predicted bounds will then allow for a unique test of the viability of thermal leptogenesis 
in the laboratory. 
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Appendix A 



Conventions and Spinor Notation 



The flat-space Lorentz metric is given by 

r/M.=??^' = diag(+l, -1,-1,-1). (A.l) 
We fix the sign of the completely antisymmetric tensor Epmpa by choosing 

£0123 = -1- (A.2) 
Greek indices z^, • • • = 0, . . . , 3 denote space-time indices. 

A.l Weyl Spinors 

A two-component complex undotted Weyl spinor (left-handed Weyl spinor) trans- 
forms in the (f,0) matrix representation of the Lorentz group S0(3, 1), i.e. under 
SL(2,C), while the dotted Weyl spinor (right-handed Weyl spinor) is in the con- 
jugate representation (O, Both spinors are related by hermitian conjugation, i.e. 
(^q)^ = and (^Q,) = • Explicitly, for a Lorentz transformation M G SL(2, C): 

C = mA/3, ~C = M*Jlp, (A.3a) 

= M-^^^C'^, e'" = (M*)"^ . (A.3b) 

^Though we make a distinction between Einstein and Lorentz indices in the beginning of Chapter 2, 
greek indices stand for flat spacetime indices in the other chapters of this thesis. 
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A. 2. Four-component Spinors 



Spinor indices are pulled by the Lorentz invariant e-tensors 



namely, 











r 





(A.5) 
(A.6) 

Furthermore, we define the Pauli sigma matrices (index 1,2,3) with lower Lorentz indices: 

The standard convention for the contraction of anticommuting Weyl spinors is 

^r? ^ = = = , (A.8a) 

1^ = = e'^^fiVa = -e'^^lpVa = ^ ■ (A.8b) 
Note the spinor index structure of the sigma matrices a . One defines 

-Maa_^a/jga/3^M (A.9) 



as well as 



^ 1 (^M^ . _ a-^^^a'^-f^) , (A.lOa) 



A. 2 Four-component Spinors 



In the Weyl basis, the Dirac 7 matrices read 
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A. 2. Four-component Spinors 



They satisfy the Clifford algebra 

{7^,7,} = 2??,,, (A.12) 
and anticommute with 75 = i7'^7^7^7^, i.e. {7/^,75} = 0. In this representation^: 

We can write a Dirac spinor in terms of a left-handed and a right-handed Weyl spinor 

= fe] , (A.14) 



V'(D) = ^(D)V = U" (A.15) 



and its adjoint spinor as 



With the chiral projectors Pl = ^{t + 75) and Pr = |(1 — 75) left-handed and right- 
handed four-spinors are given as 

and 

respectively. For the adjoints of the chiral spinors one finds tpj^ = iP(^d)Pr and 

By virtue of the charge conjugation matrix C an equivalent realization of the Clifford 
algebra is given by the transposed 7 matrices: 

C-\C = -7J, (A.17) 

with 

Ct = = = -C and = -1. (A.18) 



^The sign-convention for 75 results from the choice (A. 7) which aUowed for a clean transition in 
Lorentz signatures. 
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A. 2. Four-component Spinors 



The matrix C can be written as 



so that the charge-conjugated Dirac spinor of (A. 14) then reads 



V^^D) = Ci^l, = 1^^^ j . (A.20) 

A Majorana spinor is equal to its own charge-conjugate, i.e., tp(M) = V'{m)) that it can 
be written as 

AM)=\bA- (A.21) 
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Appendix B 

Feynman rules 



Here we provide the complete set of the Feynman rules which is necessary for the calcu- 
lations performed in Chapter 3. The method of a continuous fermion flow is addressed 
in Section 2.9. For details, see [13]. 

The momentum P always flows from the left to the right for the external lines and 
propagators shown below. Furthermore, momenta are assumed to flow towards the 
vertices. 

External Lines 




Gauge bosons A^°''^ 



Gravitinos V'/j 



65 



Propagators 



Gauginos A^") 



Matter fields xl^ 



P"^ — 

i{-f + m^) 
p2 — 



5'^ 



6'^ 



Scalar s 



P2 _ ml 



5'^ 



Gauge bosons A^°''^ 



a, II 



6, p 



p2 ^ (P2^2 



sab 



Relevant Gauge Vertices from Eq. (2.73) 




-70 T^'^^-^^^Pt 






■lV2gaT^ljPL 



iV^gaT^CiPR 



J 



^Note that for unbroken electroweak symmetry, the only non-vanishing masses for the chiral matter 
fermions arise from the higgsino mass parameter /j.; see (2.59) and (2.79). 
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i,K 



C, (7, K 




a, v, P 
o,p,Q 



Gravitino Vertices from Eq. (2.82) 




67 



Gravitino Vertices from the Effective Theory (2.83) 



Note that in the effective theory for Ught gravitinos, the external gravitinos are treated 
as Major ana fermions. 




ml - mt 




Sab \Y, 




c,p 



b, V 




igaMa 



2V6M] 




gaMa (a) 



pm; 
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Appendix C 



Hard Production Rate 



In this Appendix, we derive the hard part of the gravitino production rate for the fuU 
SM gauge group SU(3)c x SU(2)l x U(1)y- We will follow the approach of [2] but show 
in detail how to perform the integrations and point out where differences between our 
SU(3)c result and the one obtained in [2] emerge. 

Recall the definition for the production rate (3.14) 

3 



1 



hard (2vr)32i^7 4^ 



n 



3 

Y. 



/bfb I M^^i^ 1 2 + /bbf I M^'gp 1 2 + /fff I M^^)p 1 2 



{2TTf6\Pi+P2-Pz-P) 

0(|Pi -PsI 



k*) 



(C.l) 



where the shorthand notation for the products of the quantum statistical distribution 
functions is given in (3.15a). In the following, we show in detail how to obtain 



hard 



d^p 



+ 



G 



hard 



d^p 



+ 



G 



hard 



d'^p 



(BBF) 



hard 



(C.2) 



C.l Production Rate for BFB Processes 



To keep better track of the various contributions, we use the notation 



dV 



G 



d^p 



(BFB) 



hard 



1 Mp 



3m1 

G 




A/- 1 2 

''BFB, j -'bfb 



which was introduced in Eq. (3.23). The integrals I-q^b ^i^l be calculated below. The 
thermal bath provides a distinguished frame of reference. It is the frame, in which 



(C.3) 
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C.l. Production Rate for BFB Processes 



the quantum statistical distribution functions /b and /p have their simple form (3.2a). 
Hence, we perform all integrations in the rest frame of the plasma. 



C.1.1 Contribution from |-Mi| 

We have to compute the integral 
.|Mi|2 _ 1 f d^p 



^BFB - (2vr)32i5; J An 



j-j- d^Pi 



{2TTf5\Pi + P2-Pz-P) 



x/BFB|Mi|2e(|pi-p3|-F) . (C.4) 



The squared matrix element 



2^2 



|Mi|2 = -t-2s- — (C.5) 



becomes singular for t = {Pi — P^)"^ ^ 0. It proves useful to relate the phase space 
integrations to the reference three momentum k = pi — P3. The Lorentz invariant 
measures can then be written as 



2Ei 



5{Pl)Q{Ei)dEid?pi 

j dPk5^{\L + p^ - pi)S{P?)&iEi)dEid^pi 

6{Ef - |k + p3\^)eiEi)dEid^k (C.6) 



and 



^6\P, +P,-P-P,) = ^-^^^-^Q{E,)dE,d^p, 

^i^2 ^|P2| 

X b{Ex + E2-E- ^3)5^(pi + P2 - P - Pa) 

= ^(^ + ^V^^-'P + P^-P^'^e(i^ + E, - E,) 

2|p + P3 - Pll 

= 6{{E + E3-Eif-\p- k|2)e(S + ^3 - El) . (C.7) 

In order to perform the angular integrations, we are free to choose 

k = k (0,0,1)^, 

p = E (O,sin6i,cos6i)^, 

P3 = E'3 ( COS sin 0, sine?!) sin 0, COS 0)"^ , (C8) 
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C.l. Production Rate for BFB Processes 



where we have made the approximation that the typical energies in the thermal bath 
are much higher than the rest masses of the particles involved in the scattering. This 
yields 

s = [p^ + p^f = {P + P^f = 2EE^{1 - sin sine/) sin ^- cos 61 cos 6*) , (C.9a) 



t = {Pi- Pif = {El - E-i 



(C.9b) 



and allows us to write |Mi(t,s)p as \Mi{E^Ei^E2,k,0,6,(j))\^ . For the 5- functions in 

E^ + -{E + E^i-Eif' 



(C.6) and (C.7), we find 

<5((£; + S3-^i)'-|p-kp) 



(5(£;f-|k + p3| 



2kE 



2kE', 



6 cos 9 



5 ( cos ( 



2kE 

El - El - e 
2kE-x 



(C.IO) 



Equation (C.4) then reads 



BFB 2i^7r9 E'2 



dcos 

2 T?2 u2 



X 5 ( COS 6* - — — ^ — ^ 1 ( COS 



2kE'i 



dEi I dEsi-d cos e)d(j) j dkdQk 

E^ + e-{E + E3-Eif 



2kE 



X/] 



BFB 



Mi{E,EuE2,k,e,e,(^) Qik - k*)e{Ei)Q{E3)Q{E + Es - El) . 

(C.ll) 



The integrations over cos 9 and cos 9 yield restrictions on the integration range. Using 
also the 0-functions in (C.ll) and positivity of k we find 

cos9 < 1 ^ k > El - E3 , 
cos9 > -1 ^ Es- El < k < Ei + Es , 

cos9 < 1 ^ El - E3 < k < 2E + E3 - El , 
cos9> -l^k> E3-E1 . (C.12) 

The relations (C.12) are all fulfilled with the inclusion of Q{k — \Ei — E^l) and 
9(2^; + ^3 -^1 - k). 

Furthermore, it is easy to integrate out all angles and we find 



4fB = 212^6^2 / dEidE3dk f BFB 9bfb ^ 



(C.13) 



where 



Ci] - ((^1 - ^3)^ - k') 



2 EI + EI + 2EE2 {E-i + Eif{E + E2f 
3 fc2 k^ 



(C.14) 
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C.l. Production Rate for BFB Processes 



and 

n = e{k- k*)e{k - \Ei - 

X Q{Ei + E3- k)e{2E + E3- Ei-k) 

X e(^i)e(^3)e(^ + ^3 - e^) . (c.15) 

Here, E2 = E + E^ — Ei is understood. Using 

6(^1 + E3-k) = l-Q{k-Ei- E3) 

and 

e^k -El- Es) = e{k -El- Es)e{k - \Ei - Es\) , 

we can rewrite (C.15) as 

n = \e{k - k*)Q{k - \Ei - E3\)e{2E + E3- Ei-k) 
- Q{k - k*)e{k -El- E^)<d{2E + E3 - Ei - k) 
X Q{Ei)Q{E3)Q{E + E3- El) . (C.16) 

Now we insert 1 = Q{k* — Ei — E3) + Q{Ei + E^ — k*) in the second term in the brackets 
of (C.16) and use 

e(A; - k*)Q{k -El- Es)e{k* - Ei - E3) = e{k - k*)e{k* -El- E3) , 

and 

e{k - k*)e{k -El- E3)@{Ei + Es- k*) = Q{k - Ei - ^3)6(^1 + E3 - k*) . 
This leaves us with a sum of three integrals which we have to compute: 

-^BFB = 5bFB, 1 + 5'bFB, 2 + S'bFB, 3 ' (C-17) 

namely, 

SbfI 1 = ^12^ dEi dEs fBFBQ{E + E3- El) 

X J dk@ik-k*)Q{k-\Ei- E3\)e{2E + E3- Ei-k)g^^p'i\ (C.18) 

#B,2 = -^12^ dE3 fBFB@iE + E3 - Ei)Q{Ei + E3 - k*) 

X / dk Q{k -El- E3)e{2E + Es-Ei- fc)^!.™' ' (C-19) 
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C.l. Production Rate for BFB Processes 



SbfIs = -^12^ fBFB@{E + E^)Q{k* -El- E3) 

X [ dk Q{k - k*)Q{2E + E3-E1- k)g^^i'^^ , (C.20) 



Calculation of Sf^g ^ 

I A/ \^ 

The 0-functions for in the integral S'gpB i can be further manipulated. We multiply 
(C.18) by 

1 = e{k* - \Ei - Esl) + @{\Ei - Esl - k*) , 
which allows us to split the integral into two parts 

ffBFB, 1 = A'bfb, 11 + 5'bfb, 12 • 

(C.21) 



Note that 



and 



Q{k - k*)Q{k - \Ei - EsDOik* - \Ei - Esl 
= Q{k-k*)e{k* - \Ei-E3\) , 



Qik - k*)eik - \Ei - EsDQilEi - E-sl - k*) 
= e{k - \Ei - E3\)e{\Ei - Es\ - k*) , 



so that we find 



2 Q /"^ 

aBFB, 11 = ^12^ dEifFBF@{E + E; - i?i)e(A:* - \Ei - E^l) 

X J dk @{k - k*)e{2E + E3-E1- k)g^^p'^'' , (C.22) 

and 

9bfb, 12 = ^12^ dE-s dEifFBF@{E + E3 - Ei)Q{\Ei - E^l - k*) 

X j dk e{k - \Ei - Es\)e{2E + E3-E1- k)g^^i,i' ■ (C-23) 

The integration over dk in n yields a lengthy expression. From Q{k* — \Ei — E-^l) 
in (C.22) we see that the Ei integration only contributes for 

Es-k* < Ei< Es + k*. 
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Hence, in the limit k* — > 0, we can set Ei = E3 in fpBF and find 

12 1 /"OO 

9BFin = ^ME)J^ dE^Elh{E^){l + UE^)). (C.24) 

This integral can be evaluated analytically by choosing the proper series expansion for the 
exponentials which reduces the integral to a Laplace transformation on the summands. 
The transformed series can then be resummed [46]. For (C.24) we find 

yBFB.ll ~ /IQ^6 / "-^^ 



487r6 7o (e^ - 1) 



2 



T^h{E) 



^ poo 

n I dxx^e 



2 —nx 



487r6 

n=l 

T^h{E) ^ 1 T^h{E)a2) 

n=l 

where C,{x) is the Riemann Zeta function. 

Now consider ffgpg 12- ^''^ order to remove the absolute value in the B-function of (C.23) 
we insert 1 = Q{Ei — E^) + Q{E^ — Ei) and thereby split 12 parts, 

|Mi|2 _ |Afi|2 |Afi|2 
yRFB, 12 ~ yRFB, 121 ^BFB, 122 ' (,*^.ZDJ 

where 

5bFB, 121 = ^12^ dEifBFBQiE + E3 - E^)Q{E^ -E-s-k*) 

X Jdk e{k -Ei + Es)e{2E + E^-Ei- k)g\^^'^ , (C.27) 

5iBFB,122 = ^T2^y^ ^^3 di?l/BFBe(i?3 " i?l " A;*) 

X j dk e{k -E3 + Ei)Q{2E + E3 - Ei - fe)^!^^'/ . (C.28) 
Performing the k integration yields 

to, 121- dEsJ^^^^^dE.h^^^—^, (C.29) 

5bFB,122 - "28^6^2 dEi/BFB— • (<--30) 



Since 



, _ /^E^X : _ (C.31) 
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an integration by parts will extract the logarithmic dependence on k* in (C.29) and 
(C.30). Let us consider first equation (C.29) which is then given by the sum of the 
surface term and the remaining integral, namely, 



Afi 1 2 _ I Ml p surface | Mi | ^ partial 
ffBFB, 121 ~ 5bFB, 121 + 5bFB, 121 



(C.32) 



The surface term is 



|Mi p surface 
5bFB, 121 



1 



287r6^2 
fF{E) 



In 



\E'\ — E-i 



-1 Ei=E3+E 



/bfb(-£^i + E'^)E2 



Ei=E3+k* 



OO / T,* 

dEslni^ 



(C.33) 



where we used E2 = E + — Ei and set /c* ^ in the distribution functions after 
evaluating the borders. For the analytic integration of (C.33) we write the logarithm as 
\x].{k* /E^) = In A;* — \n.E^. The integration of the first term is then the same as in (C.25). 

For the second part we have to evaluate an integral of the form 
dx In X 



X _ 1)2 

00 



E 

n=l 

00 



dx In X e 



n=l 



3 27 2 Inn 



Thus, for (C.33) we find 



I Ml I ^ surface 
9bFB, 121 - 3 . 277^4 



In 



77r 



2C'(2) 



3 C'(2) 

-I ry < 2_LZ 

2 ^ C(2) 



(C.34) 



(C.35) 



where we used ^(2) = vr^/G. In the remaining integral, we can set /c* — > and obtain 



Afi p partial 

9bfb, 121 



1 



dE^ 



E+E-i 



dEi In 



\Ei — -E3 1 



Ei 



xGiE,-Es)—[fBFBEl{Ef + El)] . 
Now we turn to the second integral (C.30). Note that for the surface term 



|Mip surface 

5bfb, 122 



28 



dE-i 



In 



\Ei — E3 
E3 



fBFB{Ef + -Ef 



Ei=E3-k* 



Ei=0 



we get an additional contribution from the lower integral border. 



lim 

Ei^O 



In 



1^1-^31 

E-s 



fB{Ei){El + El) 



-E^ T . 



(C.36) 



(C.37) 



(C.38) 
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In fact, the upper border yields the same contribution as in (C.35), so that we can write 

|Mip surface | Mi ^ surface |Miprest „q\ 

5bFB, 122 ~ 5bFB, 121 + S'bFB, 122 (^U.^yj 

with (in the hmit k* — > 0) 

rj-i poo 

■ / dE^ E-s [fBiE-s) + ME + E3)] . (C.40) 
Jo 



|Mip rest 

5bfb, 122 



2%6 
Th{E) 

./o 



This can be written as with E^ = xT and a = E /T as 

IMiprest _ T^hjE) ^ r . r -nx 
5bFB, 122 — 287r6 2^ j ^ 

T^h{E) 



28 



n=l 

00 

E 

n=l 



1 



(C.41) 



and be resummed to 



9bFB,122 " 



28^6 



Lio 



The dilogarithm Li2(x) is defined as 

U2{x) -■ 



^ , hi 1-t 



(C.42) 



(C.43) 



The authors of [2] find a contribution hke (C.42) for the BFB-part of the gravitino 
production rate, but do not specify from where it emerges. In contrast, we find below 
in (C.49) that (7^3 ^22* exactly cancels out and thus we do not get (C.42) as part of the 
production rate. 



For the remaining integral, we find for A;* ^ after an integration by parts 

1 



|Mip partial _ 
^BFB, 122 - 287^6^2 



Jo 

d 



E3 



e(i?3 - E,)— [fBFBE\El + El)] . 



(C.44) 



Calculation of Sf^^ 2 

In equation (C.19), the 0-functions allow for a contribution of the integrand only if 

Ei + E^KlE + Es-Ei-^EiKE . 
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We can take the limit k* ^ and perform the k integration which gives 

9b?b,2 = dE.fBFBiEi - E) [{E + E^)E^ + E,{E - Ei)] . (C.45) 
This can be rewritten with E2 = E + E^ — Ei as 

& = -^8^^ f dE, Ij di^i/BFB^^ {EM - E'El) . (C.46) 
An integration by parts yields a surface term and we write 

„\Mi\^ _ |Afi|2 surface I Afi 1 2 partial (r A7^ 

yRFB, 2 ~ yRFB, 2 i/BFB, 2 ■ \^-^') 

For the lower border i?i = of the surface term 

El — -E3I 



^ -1 POO 

JMil^ surface _ / .p 

5bfb,2 - 287r6^2 



/bfb In 

X {E-t{E + Es-Ei)-'-E^El) 



Ei=E 

'2/771 I 771 z:^ \2 Z7i2 7:i2\ 



(C.48) 



Ei=0 

the Bose-distribution fsiEl) becomes singular. In the limit Ei — > 0, we obtain 

p rp poo 

surface ^ _^ dE^E^E + E,)[l + ^E,)] 

= -&t^^ ■ (C.49) 

We see that the arising dilogarithm (and (^-function) from the integral (C.40) cancels 
with (C.49). 

In the partial term, we include Q{E — Ei) in order to change the upper integral border 
and find 

5bFB,2'''''*'''' = o8„6 ir2 / dE^ I dEi In 



X QiE - Ei)^^ [fBFBiEfEi - E^eD] . (C.50) 

Calculation of Sf^g 3 

The k integration in (C.20) can be carried out directly and leads a lenghty result. From 
the G-function e{k* - Ei - E^) in (C.20) it follows that Ei < k* and E2 < k*. Thus, 
one finds that the integration over k yields an expression of order k* and hence 

5Bra3 = for k*^0. (C.51) 
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Result for L 



BFB 



Collecting the results, we find for the contributions from jiVfip 



'BFB, 



jMi p surface 



'BFB —•6'bFB,11+ ^^BFB,121 +5bFB,121 + 5bFB, 122 + ^BFB, 2 



+ 



1927r4 

1 



Mip partial 
1,121 

C'(2 



Ml I partial 



Ml p partial 



, , r\ 17 



rE+Ey, /If F I 



G(£;i - ^3) 



d 
d 



■^2 ^772 



/bfb-^(^i + £"3; 



+ Q{E:, - E,)— [fBFBiEl + El)] 



e{E-Ei) 



dEi 



/bfb( 



ElEl 
E^ 



Ei 



(C.52) 



This concludes our analysis of | Mi | ^ and we can turn to the integrals containing | M2 \ ^ . 



C.l. 2 Contribution from IM2I' 



The expression to evaluate is 



t\M2\^ 

'bfb 



1 



dQr 



{2'kY2E J An 



d^Pi 



n {2nf2E, 



{2n)H\Pi + P2-P-Pz) 



x/BFB|M2|2e(|pi-p3|-F) . 

Since there is no singular behavior of IM2I 



(C.53) 



: t for t ^ 0, we set /c* — > from the very 
beginning. We use the center of mass momentum q = p + ps as a reference momentum 
which yields in an analogous manner to (C.6) and (C.7) the integration measures 



2Ei 



and 



5{{E + E3- E2f - |q - P2|^)e(^ + ^3 - E2] 



2E3 



6{Ei-\ci-p\')Q{Es)dE3d-'q. 



(C.54) 
(C.55) 



Rotational invariance allows us to choose a frame where 



q = g (0,0,1) 

p = ^ (0, sin 6i, cos (9)^ , 
P2 = E2 (cos (j) sin 9, sin sin 9, cos 9)^ 



(C.56) 
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Again, we have assumed that the typical energies in the plasma are much higher as 
compared to the rest masses of the external particles which are involved in the scattering. 
The squared matrix element IM2P = t in system (C.56) reads 

= -2£;£;2 (1 -cos6lcos6i-sin6lsin6isin<?i) . (C.57) 

For the (5-functions we find 

5{{E + E,- E,f - |q - = -^5 f cos 6 + - {E + - E,f 



2qE2 V 

1 / ~ E^-El 

— -5 cos 6* ^ 

2qE V 2qE 



Equation (C.53) thus reads 



-^BFB = j {—dcos6)d(l>{—dcos6)d4>dE2dE?jdqdQ.q 

( E'+q^-{E + Es-E2r \J , E'-El + q' 
X cos 6 cos 



2qE2 J V 2^£; 

X /BFB|M2|'e(^2)e(^3)e(^ + E3- E2) . (C.59) 

From the integration over the 5-functions, we find the phase-space restrictions 

cos9 < I 2E2 - Es - E < q < E + Es , 
cose > -1 <^ q > E + Es - 2E2 , 

cosO < 1 <^ E - E3 < q < E + E3 , 

cos9> -l^q> E3- E , (C.60) 

which yield the G-functions Q{q -\E- E^l), e{E + E3 - q) and Q{q - \2E2 - E^-E]). 
We can perform the remaining angular integrations and find for Eq. (C.59) 

4fb' = ^TT^ / dE2dEsdqhFB g^^H ^ (C.61) 



with 



S = {iEs + Ef-q') 



-1 + 



E^ — 2E2E3 — E'^ + 2E2E 



(C.62) 



and 



n = eiq- \2E2 -E3- E\)Q{q - \E - E^l) 

X + ^3 - q)Q{E + E3- ^2)6(^2)6(^3) • (C.63) 
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Now we use 

e{q - \2E2 -E3-E\) = l- e{\2E2 - E3 - E\ - q) 

and 

9(12^2 -Es-E\- q)e{E + E3-q) = 9(12^2 - E3 - E\ - q) 
to split (C.61) into two parts, namely, 

-^BFB ~ 9'bfb, 1 + S'bfb, 2 ' (C.64) 

with 

X j dqe{q -\E- Es\)e{E + E^ - q)fBYB&t , (C.65) 

CB,2 = -^TT^y^ dE2j^ dE,e{E + E,-E2) 

dqQ{q -\E- E3\)Q{\2E2 - E3 - E\ - q)fBFB9^BFB ■ (C-66) 
We can take away the absolute value in Q{q — \E — E^l) if we insert 

l = e{E- Es) + 9(^3 - E) , (C.67) 
which splits again the integrals (C.65) and (C.66) into 

5bfb, 1 ~ 5bfb, 11 + 5bfb, 12 ' (C.68) 

S'bFB,2 ~ 5bFB,21 + fi'BFB,22 • (C.69) 



X 



Let us consider first (C.68) which then reads 

foo rE+Es rE+Es 



&i 11 = ^nTT^ / dE^ / ^^2/bfb / dqQ{E - E3)g^^^i , (C.70) 
9BFB,i2 = inrTW2 dEs dE^fBFB dq Q {E^ - E) 9^^^^ . (C.71) 

I^'^TT^'E^ Jq Jq JEs-E 



Integration over q gives 

fOo i-E+E-i 



= "^^""Iq^ ' "^^'^^^ " ^3)/bfb^| - , (C.72) 

Uf 12 1 f-E+Es / jp \ 

9bfb, 12 = ^8^^ I dEs dE^eiEs - E)f^^^E^ - E2J . (C.73) 
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Now we turn to ggpg 2- '^^^ insertion (C.67) yields 

5bFB,21 = ~2ll7r6^2 / ^^3/ dE2fBFB 



JO 



X y dqe{E - Es)e{q -E + E:,)Q(\2E2 - E3 - E\ - q)g^^^,f , (C.74) 

5'BFB,22 = ~2ll7r6^2 /) dE2fBFB 

X / dgG(i?3 - i^)e(g + - i^3)e(|2£;2 - ^3 - e| - g)5i^£f . (c.75) 



^BFB 

We proceed in the same manner as in (C.67) and insert 

1 = 6(^3 +E-2E2-q) + e{2E2 - E3 - E - q) , (C.76) 
which leads to four integrals 

yBFB, 21 ~ WBFB, 211 "t WBFB, 212 ' \^-' ' ) 

5bFB,22 ~ 5'bFB,221 + 5bFB, 222 " (C.78) 



Note that 



Q{q-E + E3)e{E + E3- 2E2 -q)^ E2<E3 , 
Q{q-E + ^3)0(2^2 -E-E3-q)^E <E2, 
Q{q + E- E3)e{E + E3-2E2-q)^ E2< E , 
Q{q + E- ^3)0(2^2 - E- E3-q)^ E3< E2 , 

so that we include the corresponding G-functions, and then integrate over q, namely, 
,,,,2 1 r'^ fE+Ea I-E+E3-2E2 2 

Cb, 211 = - inrYW2 / '^^3 / dE2 /bfb / dqe{E3-E2)e{E- E3)g^iii 

= 3^8^ I ^^3 dE2@{E3 - E2)Q{E - E3) 

X /bfb(^2 - E3) [{E2 - E3){E2 + 2E3) - 3{E2 + E3)E] , (C.79) 
,,,,2 1 rE+Ea i-2E2-E-Ea 2 

Cb,212 = -^TT^ dE3 dE2 /bfb J^_^ dqe{E2 - E)Q{E - E3)gZi 

1 1 j-oo rE+Ea 

= -3^8^ J ^^3 J dE2Q{E2 - E)&{E - E3)fBFB{E2 - Ef , 

(C.80) 



81 



C.l. Production Rate for BFB Processes 



ffBFB,22i = -^TT^ dE, dE2 /bfb _^ dqQ{E - E2)eiE3 - E)g\^^i 

1 1 roo rE+E:i 

= 3^8^ J J dE2Q{E - E2)e{Es - E)fBFB{E2 - Ef , 

(C.81) 

5l3ra,222 = -^TT^ / dE^I dE^fBFBi _ dqQiE2-Es)eiE,-E)g^i^i 



1 1 



^0 JEz-E 

/ d£;3 / dE2Q[E2 - E-i)e{E^ - E) 
Jo Jo 



3 28^6^2 

X /bfb(^2 - ^3) [(^2 - E3){E2 + 2E3) - 3(^2 + Es)E] , (C.82) 



Result for 

For the contributions from IM2P we thus find 

j\M,\^ _ \Ah\' |Af2p IM2P IM2P IM2P |M2| = 

-"BFE ~ i/BFB, 11^ i/BFE, 12 i/BFB, 211^ i/BFB, 212 I/BFB, 221 ^ i/BFB, 



222 



1 



2567r6 7o 



d£^3 / dE2fBFB 







X i + - ^3)|l f ^ - El) + e(^3 - ^) f f - E2 



E"^ \ 3 J ' ' V 3 

E2 — E3, 



+ e{E3 


-E2)Q{E- 




-e{E2 


-E)e{E- 


Es)^- 


+ e{E- 


- E2)e{Es - 




-e{E2 


- Es)e{E3 


-E) 



[{E2 - E^){E2 + 2E^) - 3(^2 + ^3)^] 



3£;2 

E2 - Ef 

3^2 

E2 - Ef 

3^2 

^3^2^' [(^2 - E^){E2 + 2Es) - 3{E2 + Es)E] | . 

(C.83) 

This concludes the calculation for the BFB-processes since there is no contribution from 
the matrix element |M3|2 because of cj^g 3 = 0. 
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C.2 Production Rate for FFF Processes 



We now turn to the processes where all in- and outgoing particles for the 2 — > 2 scatter- 
ings are fermions. We have to compute 



FFF,i -"FFF 



(C.84) 



C.2.1 Contribution from |Mip 

Let us first discuss matrix element |Mip. The solution of the integral 



r|A/i| = 

'fff 



{2tt)^2E J An 



{27r)^6\Pi + P2-P-P3) 



x/fff|Mi|2G(|pi-P3|-F] 



(C.85) 



is obtained as in section C.1.1. Instead of /bfBi we now have to integrate over the 
statistical factor /fff (3.15a). 

The analogous expression to (C.24) is 

f oo 



^FF, n = ^M^) dEsElMEs){l - MEs)) 
and the E-^ integration yields 



(C.86) 



T^fFiE) 



5fff,ii- 48^6 ''^(e- + l)2 



487r6 ' 



n 



n=l 



48^6 ^> 



n=l 



487r6 



(C.87) 



Now consider g^-p^ 12 which we have written in section C.1.1 as two parts 



ypFF, 12 ~ ypFF, 121 "T i/FFF, 122 



namely, [cf. (C.29), (C.30)], 



|Mi|2 _ 1 , ,p 

yFFF, 121 - 28vr6£'2 ' 



/O 



1 



|Mi|2 _ - / ,p 

yFFF, 122 — 287r6S2 ■ 



'3 / "-C/1/FFF ^ , 

E-i+k* ^1 - ^3 

rEs-k* Tp2(Tp2 i 1^2 1 



JO 



dEif, 



FFF- 



E\El + El 
El — £'3 



(C.88) 

(C.89) 
(C.90) 
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An integration by parts will exhibit the logarithmic k* dependence in the surface term 
so that we find for (C.89) 

A/ip _ lAfipsurface |Mippartial Q1\ 
9fFF, 121 ~ 9fFF, 121 +5fFF, 121 ' (u.yij 

where the surface term reads 



^ siirfarp 1 



„|A/ir surface _ . ,p 

yFFF,121 - 98^6 p2 / 



287r6^2 



oo 







E-^=Ez+E 



Ei=Ez+k* 



277r6 ^ \EsJ (e^3/T + i)2 • 

Again, E2 = E + E3 — Ei is understood and we have set /c* — > in the distribution 
functions after evaluating the borders. The integral is solved in an analogous manner as 
in (C.33) so that we find 



All I ^ surface 

9fff, 121 - 3 . 287^4 



\ [2T\ 3 C'(2) 



(C.93) 



\k* J 2 ' C(2) 

Note that in contrast to (C.35) there is a factor of two in the nominator of the logarithm. 
For the remaining part, we find for k* ^ 

IMippartial 1 ["^jip C^^^^^tj^ i ( \El - E-^ 

5fff,121 - - 08^6 Z72 / "^3 / aEi in 



X Q{E^ - i?3)^ [hYYEl{El + El)] . {CM) 

The second integral (C.90) does not yield any additional contribution from the lower 
border as in (C.38). Thus, 

= , (C.95) 



and we easily obtain 



|Mip surface |Mip surface c\p.\ 

5fff, 122 ~9'fff,i2i • (u.yoj 



Furthermore, we find for fc* — > 



|Mi Impartial 1 [°° ^rp f^^^^^ip i f 1^1 - E^ 
5fFF, 122 - o8„6 F2 / '^^3 / dEi hi 



dEi 

An integration by parts of 

foo rE 



2^7r^E^Jo 'Jo ' \ Es 

X e{Es - E,)-^ [hFFE\El + El)] . (C.97) 



&2 = I (^1^2 - E'El) (C.98) 
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does not yield surface term; see (C.49). Hence, we obtain 

1 



ypFF, 2 



oo r-E+Es { \Ei — E'A 

2^TT^E^ V ^3 

Q{E - [hMElEl - E^El)] . 



(C.99) 



LA/ P 

Again, g-p^Y 3 ~ 0, since the argument preceding (C.51) is independent of the quantum 
statistical distribution functions. 



Result for /i 



FFF 



For the FFF processes, we find in total from the contribution of 

FF 

T^h{E) 



^|Mi|2 _ |Mip „ |A/i|2 surface |A/i p partial |Mi |2 partial |A/i|2 
-'FFF ~ ^FFF, 11 + ^ ^FFF,121 + ^FFF, 121 +5fFF,122 + 5bFB, 2 



3847r4 
1 



' (2T\ 17 C'(2) 



2567r6./o 'io V ^3 







6(^1 - ^3) 



d 
d 



fBFB^{Ei + E^) 



+ eiE3 - E^)— [fBYB{El + El)] 



+ e{E-Ei) 



d 
IE 



f ( ^1^2 p2^ 
JBFB(— p 



(C.lOO) 



Note that IM2P does not contribute to the FFF processes because c^gg = so that we 
can immediately turn to the integrals containing jMsl^. 



C.2.2 Contribution from \Ms\'^ 

We now encounter for the first time also the matrix element IM3P, namely. 



|M3|2 

FFF 



d^r 



(27r)32^ J Att 



n 



d^Pi 



\ {2TTf2E, 



{2^f5\Pi + P2-Pz-P) 



x/FFF|M3|2e(|pi-p3|-F) 



(C.lOl) 



Recall from (3.21) that |Af3p = t^/s so that we can set fc* — > and employ the for- 
malism developed in section C.1.2. We choose (C.56) as the frame of reference for the 
integrations. 
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It follows that 



s = {E + E^f - , 

and with t as in (C.57) we find for the squared matrix element 
1 



2EE2 I 1 — cos 9 cos 9 — sin 9 sin 9 sin < 



{E + EsY-q^ 

The angular integrations are straightforward to compute and we get 



^FTF = 213^6^2 / dE2dE3dqf FFF 5fff ^ 



A/3 1 2 



(C.102) 



(C.103) 



(C.104) 



with 



&f={iE3+Ef-q' 



2 (2Ei - 6E2E3 + 3^3 + 2^2^ - E^ 



{E3, — E)^{Ei — E2) 



2^ 



(C.105) 



and $7 as in (C.63). The manipulations of the G-functions are the same as for so 
we can merely present the results. Analog to (C.70), (C.71), (C.79), (C.80), (C.81) and 
(C.82), we find 



5fff' 11 



ypFF, 12 



1 



5fFF,211 



5'BFB,222 



287r6^2 

1 

28716^2 
1 



dE-, 



dE^ 







E+E3 ^2^2 
dE2Q{E - £^3)/fff-^— 

E+E3 E'^E'^ 
dE2Q{E3 — -E)/fff- 







E+E3 



E + E3 
dE2Q{E - E3)e{E3. - E2) 



287r6S2 

X /fff(^2 - ^3) [^2(^3 -E)- EsiEs + E)] , 
1 



dE^ 



E+Ez 



dE2Q{E3 - E)Q{E2 - E3 



while 



287r6^2 

X /fff(^2 - ^3) [MEs -E)- EsiEs + E)] 



5fff, 212 ~ , 
5fff' 221 = . 



(C.106) 
(C.107) 

(C.108) 
(C.109) 

(C.llO) 
(C.lll) 
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Result for /i 



FFF 



The contribution of matrix element IM3P yields 
^FFF ~ ypFF, 11 ypFF, 12 "i" ^FFF, 211 "i" yRFB, 222 



1 



2567r6 



00 rE+Es 

dE^ / (ii?2 /fff 

JO 



2p2 



1 i^fi?. 



E + E3 



+ 9(^3 - ^) 



1 g^^l 

E'^ E + E3 



e{E - E3)e{E3 - ^2) ^' [^2(^3 - ^) - ^3(^3 + E)] 



+ 6(^3 - ^)e(^2 - ^3) 



^2 

-£'2 — £-3 
^2 



[^2(^3 - ^) - ^3(^3 + E)]} . (C.112) 



C.3 Production Rate for BBF Processes 



The BFB processes contain the contributions with two incoming bosons and the gravitino 
and another fermion as outgoing particles. In total, 



Note that all c- 



(a) 

BBF, 1 



(BBF) 
hard 



a=l P \ G/ i=l 



I ''BBF, i -"BBF 



(C.113) 



so that we are left with the computation of IM2P and IM3 



C.3.1 Contribution from [Map 

We have discussed in great detail how to perfom the integrations for the non-singular 
matrix elements IM2P and |-/Vf3|2 in sections C.1.2 and C.2.2, respectively. Again, all 
that changes is a different combination of quantum mechanical distribution functions. 



r|M2| = 

'bbf 



d^lr, 



[lirflE J 4tt 



n 



d^Pi 



\ {2TTf2E, 



{2t:)H\Pi + P2-P3- P)/bbf|M2| 



(C.114) 
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4bf = T^^r-Q I dEs I dE2fBBF 



so that we can immediately present the result: 

r-OD fE+Ez 

2567r6 7o JO 
x\q{E- i?3)|| - E^ + e{E, - i?) (f - ^2 

+ 6(^3 - E2)Q{E - ^3) ^3~2^' [(^2 - ^3)(^2 + 2E^) - 3(^2 + ^3)^] 

-Q{E2-E)Q{E-E,) ^^'~J^' 

{E2 - Ef 



+ e{E - E2)e{E3 - E) 



- Q{E2 - ^3)6(^3 - [(^2 - ^3) (^2 + 2^3) - 3(^2 + ^3)^] I • 

(C.115) 

C.3.2 Contribution from iMgp 

For the squared matrix element IM3P we can set as usual k* — > 0, i.e., 
IM3P _ 1 [ dOp 



,=1 (27r)32i?, 



(27r)454(Pi + P2 - P3 - ^')/bbf|M3|2 



(C.116) 

and find in complete analogy to section C.2.2: 



^^BF' = ^y^ dEsJ^ dE2hBF 

( 1 s?^! 



00 rE+Es 
dEs 

Jo 



x{e{E- E3) 



E'^ E + E-i 



- e(^ - ^3)6(^3 - E2)^^^^ [^2(^3 -E)- EsiEs + E)] 

+ e(E3 - E)e(E. - E,)^ IE,(E, - E) - E,iE, + eA . (0.117) 
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Published Work from this Thesis 



The new SU(3)c x SU(2)l x U(1)y result for the colhsion term and its imphcation on 
gravitino dark matter scenarios allow for a collider test probing the viability of thermal 
leptogenesis. As a summary of this thesis, the proposed method has been published in 
Physical Review D [47] . 

An e-print of the paper is available on the arXiv server: 
http : //lanl . arxiv . org/abs/hep-ph/0608344 
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